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ABSTRACT 


The Fermi Surfaces in Cu, Ag and Au are now known to be greatly distorted,. 
with thick ‘necks’ passing through the zone boundaries. In this paper we: 
enquire whether such an electronic structure is quantitatively consistent 
with the observed transport coefficients. The mathematical model is quite: 
simple; the shape of the Fermi surface is made to depend on a single parameter 
which can be interpreted as the pseudo-potential of the {111} atomic planes 
acting on an orthogonalized plane wave, giving rise to an energy gap of 
5-10 ev at the zone boundaries. Various integrals over the Fermi surface 
can then be evaluated by elementary methods, and compared with the 
corresponding experimental quantities. The electronic specific heat and 
optical mass in the pure metals are consistent with the model. The galvano- 
magnetic effects are shown to depend a great deal on the anisotropy of the 
electron relaxation time, whose variation with energy is also probably the 
main determinant of the sign of the thermoelectric power. A better theory 
of electron—phonon interaction is needed before this, and the electrical and 
thermal conductivities, can be calculated accurately. Generally speaking 
there is no evidence which directly contradicts the rigid band model, except 
perhaps the effect of alloying on the optical absorption edges and on the 
electronic specific heat, but there are still many experimental and theoretical 


gaps in our knowledge. 
B2. 
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$1. Tae NEED FOR A QUANTITATIVE THEORY 


Tue free-electron model of a metal has an interesting history. Once 
upon a time it threw a brilliant light upon the properties of metals, and 
was the basis of the whole theory. Then the ideas of Bloch and of 
Brillouin were developed and attempts were made to calculate everything 
in terms of zones and bands and localized atomic orbitals. Now more 
recent work has suggested (Phillips and Kleinman 1959, Harrison 1959) 
that a single orthogonalized plane wave, with nearly the energy of a free 
electron, may be a remarkably good representation of the electron wave 
function except for wave vectors very near to the boundary of a Brillouin 
zone. Thus, we may return to a Fermi surface approximating to a 
sphere, in the extended zone scheme, with modifications and changes of 
connectivity only in the neighbourhood of the zone boundary planes. 

If we apply this principle to the monovalent metals, we expect to find 
very simple Fermi surfaces; a perfect sphere with one electron per atom 
does not even touch the zone boundary in the body-centred and face- 
centred structures. This is not what is observed. Na, true, seems to 
be very close to the free-electron model, but even for the other alkali 
metals there is indirect experimental evidence (Cohen and Heine 1958, 
Garcia-Moliner 1958a, Ziman 1959a) that the Fermi surface is quite 
distorted, and perhaps touches the zone boundary, in Li, Rb and Cs. 
Unfortunately these metals are not easily studied: it is difficult to make 
single crystals, there is a phase transformation at low temperatures, and 
‘doping’ with suitable impurities is not easy. For the noble metals, 
however, the evidence is incontrovertible: the anomalous skin effect 
(Pippard 1957), the de Haas—van Alphen effect (Shoenberg 1960), high 
field magnetoresistance (Alekseevskii and Gaidukovy 1959, Priestley 
1960), cyclotron resonance (Langenberg and Moore 1959), ultrasonic 
attenuation (Gavenda and Morse 1959), all show that the Fermi surface 
touches the zone boundaries in the [111] directions, with an appreciable 
area of contact. The magnetic phenomena are, indeed, absolute in that 
they show the existence of small minimal cross-sections at the Fermi 
surface corresponding to cuts through the ‘necks’ passing through the 
zone faces. One can also demonstrate that the Fermi surface must be 
multiply-connected, for electrons in ‘open’ orbits are also observed. 
These methods are capable of giving ‘topological’ data about the Fermi 
surface, and show that is very far from a simple sphere. 

The information from all these experiments is very valuable, but it 
has some serious limitations; in every case the specimen must be 
extremely pure, and must be at extremely low temperatures in order that 
the relaxation time of the electrons, for scattering by lattice waves or 
impurities, shall not be too short. One is inevitably restricted to a small 
region in the corner of the phase diagram of the metal and its alloys. 
This is unfortunate, for we have no guarantee that the electronic energy 
surfaces in k-space remain unchanged when other metals are added to 
a noble metal, or when the lattice is strongly perturbed by thermal 
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vibrations. Indeed, Cohen and Heine (1958) have specifically proposed 
that this ‘rigid band’ hypothesis must be abandoned if we are to under- 
stand certain properties of the alloys of the noble metals. It would 
obviously be valuable to have an independent check on this point. 

At the same time, there are many published measurements of the more 
‘ordinary’ transport coefficients of the noble metals. Thus, we have 
data on the electrical and thermal conductivities, thermoelectric power, 
Hall coefficient, magnetoresistance, residual resistance per impurity 
atom, etc., at all temperatures and for a wide variety of specimens—pure 
metals, single crystals, cold-worked wires, alloys, etc. There are also 
observations of the optical properties, the Knight shift, and of the 
electronic specific heats, which are not confined to the special conditions. 
of the ‘topological’ techniques. 

Now all these experimental results have one feature in common—they 
do not fit at all neatly into a free-electron model. The coefficients are 
often of the wrong sign, show anomalies of functional behaviour with 
temperature and impurity content, do not satisfy similarity conditions 
like Matthiessen’s Law and Kohler’s Rule and are certainly very seldom in 
quantitative agreement with theory. Of course we have now an easy 
way out; we attribute the anomalies to the complicated shape of the 
Fermi surface. Thus, we say that the positive thermoelectric power is 
due to the ‘hole’ states near the contact region, and that the large 
magnetoresistance is due to the large distortion of the Fermi surface. 
But it is very necessary to make these arguments quantitative and precise 
if we are to construct a consistent picture of the electronic states in the 
metal. Not all the phenomena. are even properly understood in 
principle, so that we do not know whether we can get by with the rigid 
band model or whether we must, indeed, invoke the Cohen and Heine 
hypothesis in order to explain what we see in alloys. The matter is still 
very much open to discussion. 

The main purpose of this review is to give a quantitative numerical 
framework to these problems. We use a particular model of the energy 
surfaces in the noble metals, soundly based on theory, consistent with 
the topological results, but sufficiently simple that many important 
physical parameters can be computed by schoolboy calculus. Each of 
the properties is reviewed in turn, and comparison made between the 
theory and the actual observations. In this way, one can bring out the 
really difficult points, and show where there are still serious discrepancies. 
It is hoped that a quantitative theory may encourage further careful 
and reliable experiments; it is not at all easy to be sure that some of the 
reported phenomena are not due merely to accidental impurities.in the 
specimens; for example it would be extremely profitable to make 
measurements of a number of different transport properties on the same 
specimen, at the same temperature, to avoid random errors from this 
source. A number of interesting theoretical problems are also raised. 
Near the ‘necks’ the electron energy surfaces are hyperboloidal, rather 
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than ellipsoidal, and this requires an entirely new technique in such 
problems as scattering by impurities. It will be suggested that much 
may depend on the anisotropy of the relaxation time of the electrons; 
we need a good theory of this effect. These problems are not treated 
in this review, but obviously merit study. 

I have one further, personal, motive. Having just written a book 
(Ziman 1960—to be referred to as E & P) on the theory of the transport 
properties of metals in general, I am interested in applying this theory 
in detail to a particular group of metals and, where necessary, modifying 
and extending the general argument. To avoid unnecessary formal 
theory and its bibliographical background, reference will be made to this 
book for many points which are there completely documented. 


§ 2. THE OBSERVED FERMI SURFACE 


The noble metals crystallize in the face-centred cubic structure with 
the lattice constants given in table 1. From the atomic volume alone, 
using the simple free electron model, one can deduce certain theoretical 
parameters, which are also recorded here for reference. 


Table 1 
Metal Cu Ag Au 
| ea Mine eae 3-608 4-078 4-070 
Atomic radius (A) 1-41 1-59 1-59 
Fermi energy, for free 7.04 5-51 5-51 


electrons (ev) 


Electronic specific heat 
parameter (y) for free z 


502 Re “ps 
| electrons (millijoules/mol oe 0-644 0-644 
(°K)?) 
Hall constant for free elec- 
trons = (Nec)~! (em3/ —7-45 x 10-5 | —10-65 x 10-5 | — 10-60 x 10-5 


coulomb) 


The reciprocal lattice is body-centred cubic, with the Brillouin Zone 
shown in fig. 1. This has eight hexagonal faces, with normals in the 
{111} directions and six square faces which are the {200} planes. For 
our purposes, the most important scale parameter in this figure is p, the 
distance ['L from the origin to the midpoint of a hexagonal face. In 
terms of this distance, 


DX = (2/\/3)p=1-155p; TK =TU=4/(3)p=1-225p; 
DW = \/(3)p=1-290p, 
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and the side of one of the squares or hexagons is V(2)p=0-817p. Thus L 
is the symmetry point on the zone boundary nearest to the origin, 
followed by X, K and U, and W, in that order. 

The volume of the zone is (32/34/3)p?=6-16p3. If there is but one 
electron per atom, the Fermi surface must enclose just half this volume. 
If the electrons were free, the sphere would have a radius 


(4/74/3)'8p = 0-902p. 
This sphere approaches, but does not touch the zone boundary even at 
the centre of a hexagonal face. Thus, it takes a substantial perturbation, 


Fig. 1 


The Brillouin zone for a face-centred cubic lattice, showing a sphere 
containing one electron per atom. 


of the single orthogonalized plane waves to change the noble metals 
from ideal free-electron systems; there is no direct geometrical necessity, 
as in a polyvalent metal, for the Fermi surface to cut the zone boundary 
at all. 
The actual situation seems to be as in fig. 2. The Fermi surfaces in 
all three metals are drawn to the zone boundaries, and make contact. 
The diameter of the ‘neck’ formed between the surface in one zone and 
its image across the zone boundary is quite appreciable; Shoenberg 
(1960) gives the values in table 2 from his direct measurements of the 
de Haas—van Alphen periods of orbits encircling these necks. These data 
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are consistent with the results of ultrasonic-attenuation experiments 
(Gavenda and Morse 1959) and with the magnetoresistance work, on 
single crystals (Alekseevskii and Gaidukov 1959, Priestley 1960). 


Fig. 2 


A sketch of the Fermi surface of a noble metal, showing * neck’ and ‘ belly ’ 
orbits. 


Table 2 


Cu Ag 


Neck radius, + 0-18p 0-13p 


Effective masses, my* 1:30-1:38 0-7 1: 
(Free-electron mass units) (belly) (belly) 0- 


2 (belly) 
5 (neck) 


Other finer details of the shape of the Fermi surface are still in dispute. 
Thus, the ‘belly’ circuits in the {111} and {103} planes, as observed by 
the de Haas—-van Alphen effect, do not seem quite to agree with Pippard’s 
measurements (1957) on the anomalous skin effect. The experiments 
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also give rough estimates of the ‘effective masses’ of the electrons on 
some of the orbits—i.e. the cyclotron masses 
tat Oe dA (k,) 
2r dé 
where .¥ is the area of the cross section, normal to the magnetic field, 
of the surface of energy &. For Cu, these dynamical data are more or 
less consistent with the direct cyclotron resonance measurements of 
Langenberg and Moore (1959). It is obvious that further experiments 
using these various techniques will soon make these data more precise, 
and perhaps give a definite shape to the Fermi surface in each of the 


noble metals; for the moment we shall accept the figures in table 2 
without comment. 


(2et) 


§ 3. THe ENeRGy SURFACES 


Our problem is to find an analytical function for surfaces having the 
shape shown in fig. 2, with the parameters of table 2. It would be 
agreeable if we could simply take over some calculation of the electronic 
structure from first principles, but (HK & P, p. 112) the attempts that 
have been made so far do not accord at all with the experimental results 
that we have just considered. We must fall back on some phenomeno- 
logical model which takes certain parameters from observation, and 
interpolates suitable surfaces between them. 

There is a well-known interpolation formula due to Slater and Koster 
(1954), which is an explicit analytic function, with the proper symmetry 
and periodicity at all points in the repeated zone scheme. Garcia-Moliner 
(1958b) has shown that the surface proposed by Pippard (1957) can be 
fitted quite closely by a few terms of this formula. However, this formula 
has disadvantages for numerical computations, and its theoretical 
interpretation, based on the tight-binding approximation, is not very 
sound (although it can be thought of simply as a Fourier expansion of 
the energy in k-space, without physical significance in detail). 

The contemporary philosophy of energy-band theory is quite different ; 
free electrons have come back into fashion. Our model is going to be 
based on the ‘nearly-free electron’ approximation, which would have 
us write 


LU ace> Cue OT et eae er n aeee. 13-1) 


where g are the reciprocal lattice vectors, and the ¢,,,(r) are basically 
plane waves. In the original formulation of this method (E & P, 
§ 2.2), they were simple plane waves, : 


Seag(") =expli(k+g).r],  . . . - - (3.2) 


and one found the coefficients «,,, by substituting in the Schrodinger 
equation, and solving the set of linear equations 


g 
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for all values of g. Here, of course, &, is the energy of the electron 
state y,, &,° is the energy of the unperturbed state ¢,,, and 


i : 
U,= 7 | Mees (ig) ar S20, area ie 


that is, it is a Fourier coefficient of the lattice potential which is causing 
the mixing of the simple plane wave states. 

Although (3.1)-(3.4) constitute a formal Fourier expansion of the wave 
functions and a corresponding Fourier transform of the Schrédinger 
equation, this was thought to be an unsatisfactory method in practice 
because the series (3.1) converges so slowly. There is obviously little 
value in an interpolation formula for 6, which depends on a large number 
of coefficients like Y, through the solution of a secular equation of high 
order. But this judgement is too pessimistic. 

The justification of the method comes from the theory of energy bands. 
It has long been known that in most metals a single orthogonalized 
plane wave is quite a good approximation to the wave function for values 
of k well inside the zone. In other words, if we take for ¢,,,(r) a 
suitable set of OPW’s, instead of simple plane waves like (3.2), the 
expansion of #, in (3.1) may converge very rapidly. Each OPW has 
already the rapid oscillations required to fit it into the atomic core with 
the correct number of nodes; there is no need to bring in a large number 
of higher Fourier components in order to represent this part of the wave 
function. We expect most of the coefficients that would then appear 
in place of (3.4) to be quite small and easily manageable by perturbation 
methods. 

This point has been tested explicitly by Phillips and Kleinman (1959) 
who have shown that the orthogonalization condition is equivalent to a 
repulsive pseudo-potential which is spherically symmetrical and which 
balances out most of the real attractive potential due to the ordinary 
Hartree-Fock field of the ion. This cancellation seems to be almost 
complete inside the core; only a few of the lowest Fourier components, 
corresponding to the long range part of %(r), are at all important. 

In our monovalent metals we know that the Fermi surface keeps well . 
away from all the zone boundaries except the {111} planes. This means 
that the effects, on the energy, of all coefficients of the lattice potential 
except %,, will be small. We can approximate to these effects by 
perturbation theory, writing 


|%.)? 


¢ e0 
(@) ~é& _ ea eee RY Neo 

k k ia 0 g = 
g4{111} Ok+g°— OK? 


(3.5) 
toughly speaking, this gives rise to an increase of the effective mass by a 
more or less constant factor, and leaves us still with a nearly spherical 
Fermi surface. The really important parameter in determining the 
shape of the surface is %,,,, which causes the energy splitting on the {111} 
faces at the zone, and hence fixes the amount of contact with those faces. 
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Thus, if we put into (3.1) a set of OPW’s in place of the simple plane 
waves (3.2), then we may solve (3.3) for our energy function with only 
U3,, and perhaps Wyo) not negligibly small. This should give us a range 
of interpolation functions, with only these two parameters to be fixed by 
appeal to experiment. 

Even this is too complicated for our present purpose; we wish to 
represent the Fermi surface by a simple analytical function depending 
only on a single parameter whose value will be fixed by the size of the all- 
important ‘necks’ passing through the {111} faces. It is obvious that 
U1, is the key to this. Moreover, in the neighbourhood of a particular 
hexagonal face—the one bisecting the reciprocal lattice vector g, say— 


Fig. 3 


The energy gap at the zone boundary, 


we may neglect all the terms in (3.1) except the two which are nearly 
degenerate whenk~4g. Thesecular equation then reduces to a quadratic, 
whose lower root lies below &,°, by the amount 


Oo y= 4VU(Cxn_s° = o4 Joe 4U ?} i aCe —é&,°). . (3.6) 


At the zone boundary itself, this tends to |%,|, so that the Fourier 
coefficient is just half the value of the ‘energy gap’ (at this point in 
k—space) between the two bands that are split by the action of the 
lattice potential (fig. 3). But as we go away from the zone boundary, 
we soon find that this expression is well approximated by the perturbation 


_ formula (3.5). 
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This formula is the basis of our model. We assume that 6, deviates 
from the spherically symmetric function &,° by the amount (3.6), where g 
is the reciprocal lattice vector corresponding to the nearest hexagonal face of 
the zone. The depression of the energy levels is thus greatest at points 
near the centres of the zone faces, so that the energy surfaces bulge out 
in these directions. In the neighbourhood of each {111} face, they are 
surfaces of revolution, with fairly simple geometrical and analytical 
properties. The shape of the necks depends only on the parameter 
U,, i.e. U;, in (3.6), which is therefore fixed, for each metal, by the 
radius 7 in table 2. Thus, our model has no adjustable parameters to 
tempt us into melting the theory and pouring it into the mould of the 
facts! Similar models have been used previously for the analysis of 
particular phenomena (Jones 1957, Ziman 1959b; Suffezynski 1960), 
but this is the first comprehensive study of all the transport properties 
on the same basis. 

There is one further mathematical approximation that is extremely 
convenient. The part of the Fermi surface that ‘belongs’ to a particular 
hexagonal face is certainly a surface of revolution, but its boundary— 
i.e. the frontier between two such regions—is a complicated curvilinear 
polygon. It would obviously be much more satisfactory if this boundary 
were simply a circle, as if each region of the Fermi surface were the cap 
of a circular cone with vertex at the centre of the zone and axis to the 
centre of a hexagonal face. Unfortunately one cannot map out the whole 
47 solid angle in this way with circular cones facing along the {111} 
directions; some parts are left uncovered whilst other parts lie in inter- 
secting cones (fig. 4). But in those regions where cones intersect the 
Fermi surface is nearly spherical, and its transport properties will not be 
very different from those of a portion of a sphere. The same is true of 
the regions which are left out (see fig. 5). Thus, so long as we make the 
total cap area and total volume of our eight cones equal to the total 
area and volume of the surface, it will cause very little error if we count 
some of these border lands twice over and leave out others altogether. 
The overall cubic symmetry of the figure saves us from any difficulty 
about the exact direction into which these bits of surface face. In other 
words we choose the vertex angle of each cone so that it subtends 
one-eighth of a whole solid angle, i.e. 7/2 steradians, and ignore the fact 
that they will not all pack properly together without overlap. The 
enormous gain in mathematical simplicity more than balances the small 
errors introduced by this approximation. 

As we have seen, there is good theoretical ground for our present 
standpoint. Nevertheless, (3.6) is only an approximation, and one 
might wish to improve the model for a very exact numerical computation. 
For example, what happens when we are near a [110] direction—the 
point A in fig. 1. We are in the neighbourhood of two hexagonal faces, 
both of which will exert their influence on the energy surfaces. For 
points in the plane of symmetry containing the zone edge KW, one can 
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Fig. 4 


(4) 
Eight cones cannot be fitted together without intersecting. 


Il 


12 J. M. Ziman on the 


solve the secular equation; it turns out that one must double the value 
of %,? in (3.6), and hence increase, by a factor 1/2, the energy gap there. 
As shown in fig. 5, this pulls the surfaces out a bit towards these edges, 
without materially changing them from a nearly spherical form. 

Again, we expect the [200] component of the equivalent potential to be 
comparable with Y,,,._ This will cause the surfaces to bulge a little towards 
the square faces of the zone—the regions which are not covered by the cones 
in fig. 4. But the distances from a Fermi sphere to these faces weakens 
this effect. As shown in fig. 5 (where we have taken Wgo9 = W111), there is 


_Energy contours in a section of the zone. 
as used in the 8-cone model; —- — — — allowing for the effects of the 
two {111} faces near K and of the {200} faces near X. 
The regions where cones interest, or that are left out, are shown on the right. 


little deviation from a spherical form, and the calculated value of a trans- 
port coefficient would not depend very much on the value of Yano. 
Incidentally, this may be the clue to the difficulty noted by Shoentert 
(1960), that the belly circuits normal to (100) seemed to have a larger area 
than those normal to (111). By increasing %go9, one could pull out the 
surfaces a bit towards the square faces, where they are traversed by the 
former circuits (see fig. 2), without altering the shape in the other regions. 


Ordinary Transport Properties of the Noble Metals 13 


Finally, we must say something about the form of &,°, the energy of a 
single OPW before it has been influenced by the pseudo-potential of the 
lattice. This is not necessarily the energy of a simple free-electron, 
although it is probably spherically symmetrical. We shall assume 

CUS aM cen sud fo saa ee (327) 
where m* is an ‘effective mass’ parameter, which is not known a priori. 
In other words, we introduce this as an adjustable parameter, which fixes 
the energy scale of the band structure in each metal. However, the 
experimental data are sufficient to determine m* in various ways, so that 
a consistency check is possible. In principle, m* can be calculated directly 
if one knows the core potential etc, but this is a lengthy and difficult 
computation except for an atom of small atomic number. We have 
already noted in (3.5) that the apparent value of m* will depend also on 
the lattice pseudo-potential, Y%,, as a perturbation effect from all the 
distant zone boundaries. We shall, of course, always quote m* in units 
of the free-electron mass. 


§ 4. Neck DIMENSIONS AND ENERGY Gap 


The differential geometry of the cones is surprisingly simple. Each 
cone must subtend the solid angle 47/8, so that its vertex angle is cos~!(?) 
(fig. 6). The axis of the cone is a vector, $g, to a [111] face, and this we 
take to be of length p, and as the local z axis. It is convenient to introduce 
dimensionless coordinates x=k,/p, y=k,/p, z=k,/p, and to measure energy 
in terms of the dimensionless variable 


e= &,|(h*p?/m*). beers Wes AED) 


For the Fourier component of the lattice pseudo-potential we introduce the 
parameter 


Maar Oop] Nin) Boe ee ae |. (4-2) 
The energy surfaces defined by (3.6) and (3.7) are now given by 
e=3(x?+y") +f (2), URSAE od Veoh encore ee 415) 
where 
f(z)= 44+4(1—-2)?— fu? +(l—2zy}. ww (44) 


Each surface terminates on the sloping sides of the cone at some 
ordinate z,. Since the cone has the equation 


a? + y? = 2? tan? (cos! 3) =42?, NUS Gate peste .0 ) 


we can find z,, for a particular value of «, by the relation 


e=fet+f(z). il Perey cedar eeeer (43:6) 
This same energy surface will either cut the z axis at some point z, given by 
e=f (22), eC ls eA Bn ee re Te EB) 


or it will meet the zone boundary, at z,=1, on a circle of radius 
p=prs/(2et+2u—1). . . . « + + - (4.8) 
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Our first problem is to calculate the volume lying inside the cone 
below a certain energy surface. But this is just the sum of the volume 
of aright circular cone, of height z, and semi-angle cos~ (¢)—.e. (777/27) ae 
—and the volume of the surface of revolution lying between z, and 2g, i.e. 


n[ @+eyde=20 | lef @}ae fe Ete ae Nal ee 


These volumes are, of course, in units of p*, and we must multiply 
this sum by 8 to get the total volume inside the energy surface «. We 
know that the whole zone, of volume (32/34/3)p*, would hold exactly two 


Fig. 6 


P 


electrons per atom, so that the number of electrons per atom inside this 
surface is given by 


34/3 7 “2s 
a 16 8.0[ spa | Pe—1—(—2)+2vfut+(L—2)hde | 


CA | 


34/3 7 
5 7 E 2° + (2e —1)(2,—2,) —4{(1 —2,)® — (1 —24)} 


+ [(1 —2,)+/{u® + (1 —21)?} — (1 — 29) /{u? + (1 —2q)}] 


+u*in {qa eal 4 10 
(La) b aus come ele ee ee 


We have given this formula in full to show that the integrations are 
quite elementary and all the algebraic operations are easily handled on a 
desk calculating machine. There is only one difficulty: z, and z, depend 
on ¢ through (4.6) and (4.7), so that we cannot just fix m and find w and e. 
It is necessary to proceed in the following order: fix w and choose a value 
of z,; calculate ¢ from (4.6); solve (4.7) for 2, (if z,< 1); substitute for w 
2, 2, and ¢ in (4.10) and hence find n. ; 

In this way we can construct a table of corresponding values of « and n 
for a given value of w; we have found the relation between the energy and 
the volume of the surfaces for a given lattice potential. In the pure 
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metal, »=1; by graphical interpolation one can find the corresponding 
value of e, and, from (4.8), the radius of the neck at the zone boundary. 
As we increase uw, we get a series of surfaces of this sort, all of the same 
volume but more and more distorted. Thus, we establish a relation 
(fig. 7) between the lattice potential, w, and the size of the necks on the 
Fermi surface of an exactly monovalent metal. 

It will be seen that the radius r is very sensitive to the value of the lattice 
potential, and increase rapidly once contact has been established. For a 
given value of w, it is also quite sensitive to the number of electrons in the 
system, as indicated by the curve for n=1-1,—it might be for an alloy of 
a noble metal with a polyvalent metal. 


Fig. 7 
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Calculated variation of ‘ neck’ radius r with energy-gap parameter wu, 
for 1 and for 1-1 electrons per atom. 


Now the de Haas—van Alphen experiments have given us quite accurate 
estimates of r in the pure metal. From table 2, and fig. 7 we can read off 
values of u, and hence, with (4.2), estimate the Fourier component W%,);. 
These data are shown in table 3. 

We notice at once that the appropriate values of u are not very different 
in the three metals and that they correspond to quite large splitting of the 
energy levels; the depression of energy at the zone face is between § 
and 2 of the kinetic energy of the electron at that point. These values 
of w can be translated into absolute electron energy by the definition 
(4.2)—but this requires a knowledge of m*. Yet even taking for m* 
the free-electron mass, we find that 2%,,,, the energy gap at the zone 
face, must be of the order of 5-7 ev in all three metals. If, as will appear 
later, m* is actually smaller than 1, these estimates must be raised further. 


P.M.S. c 
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This is a far cry from the values—all less than 1-5 ev—proposed by Cohen 
and Heine from the splitting of the s and p states of the free atoms. Our 
values are, indeed, closely in support of Pippard’s theoretical estimate of 
7-5ev required to give the shape of Fermi surface which he observed in 
Cu. The experimental evidence from optical absorption is not at all 
clear. Biondi and Rayne (1959) talk of a gap of 2-6ev in a 30% CuZn 
alloy and their data show this moving up to nearly 5ev in pure Cu. 
Suffezynski (1960) has made a more detailed calculation, allowing for the 
density of states in the lower and upper parts of the conduction band, 
but finds difficulty in a comparison with experiment. 

In §2 we noted that some of the ‘topological’ experiments give 
information about the cyclotron mass on various orbits. Values 
calculated from our model can be compared with some of these data, if 
we write, instead of (2.1) 


h (dk 
LAs foie a . . . . . a ee (4.11) 
Table 3 


Cu Ag 


0-180 0-126 


= 


0-20 0-168 


WM (ev: m= —= 1) 3°5 2-3 


where v, is the electron velocity normal to the magnetic field. Now 
the components at the electron velocity are, from (4.3), 


(0 20) = 7 (2. Ze), laa sa Ct) 


; m* Oz 
ie. 2 is replaced by f’(z). 
It is obvious that a ‘neck’ orbit, in the plane z=1, gives very simply 
mM,*=m* ; rye mers brea, co i Se 


this experiment should determine the true ‘effective mass’ parameter for 
our system. A ‘belly’ orbit, however, in a [111] or [100] plane, passes 
near, but not through, the axis of several of the cones (see fig. 2). Thus 
v, will have strong contributions from v,, for which (4.11) and (4,12) 
suggest 

My em metlT (als. wag) sonata are eae PT 


This function which is plotted in fig. 8 (for w=0-2) rises steeply near 
z=1, i.e. at the zone boundary. Referring to fig. 5, we see that these 
orbits cut across the zones in the range 0:8<z<0-9, so that the ratio 
in (4.14) could lie between 1-5 and 2-5. This is in rough accord with 
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the quoted ratio of belly to neck masses in Cu and Au (i.e. about 2). But for 
@ precise comparison between theory and experiment we must use a more 
complete model (e.g. including Y,,) and make more exact calculations. 


Fig. 8 


§ 5. DENSITY OF STATES 


Several properties of the bulk metal depend mainly on the density of 
states in energy at the Fermi surface. This is given in general (KE & P 


§2.11) by the integral 
N(€)= 5 | = Sones be oO.) 


()) 


where v is the electron group velocity on the element d of Fermi surface. 
For our model surface we have (4.12), which a little elementary geometry 


WE) = ae | & eM Ne SMe (82) 


C2 


gives 


18 J. M. Ziman on the 


Fig. 9 
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Density of states, as a function of electron density, calculated for various 


values of wu. The peak occurs at the point where the Fermi surface 
touches the zone boundary. 
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Density of states, as a function of wu, for 1 electron per atom. 
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For comparison we refer to the value of the density of states for a system 
of electrons of mass m* on a Fermi sphere of radius &,; 
é 1 m*k 
Mle) =e Se ey secs iin (0.5) 
Thus, the ratio of (5.2) to (5.3) can be calculated very easily from the 
values of z, and z, already obtained in the previous section, and from 
%,=0-902p. This is the quantity y/y, plotted in fig. 9, as a function of n, 
the number of electrons per atom, for various values of the parameter uw, 
and again in fig. 10 as a function of w, for n=1. Because v vanishes 
at the zone boundary there is a singularity in the integrand of (5.1), but 
this can be integrated, and all that remains is a discontinuity of slope, 
with a vertical tangent, as a function of n. 

We have used the symbol y because the most direct measure of the 
density of states is the electronic specific heat. This is not an easy 
quantity to measure accurately, but the best data for the purest specimens 
are set out in table 4 (Corak, et al. 1955, Rayne 1957a, 1958, Guthrie 
1959). This observed value is compared with the value y,, calculated 
for a simple free electron model with m*=1; this ratio is often called 
the thermal mass. 


Table 4 


Yops (millijoules/mol (°K)?) 


yr (calculated for free electrons) 0-502 0-644 0-644 


M,=Yobs/VF 1:37 0-95 1:15 

Yeale/Vo from fig. 9 \ 1-57 1-54 1-56 

Hence 0-87 0-62 0:73 
: YF Ycale 


Now from table 3 we can look up values of wu for the three metals, and 
from fig. 10 extract corresponding values of the ratio y/yo. The two rows 
do not agree very well. However, there is still the unknown parameter 
m*, which occurs in (5.2)—i.e. yo/yp=m*. We can adjust our 
theoretical values of y and make them agree with experiment by giving 
m* the values shown in the final row of table 4 (which we shall label 
m*,). These values are more or less consistent with the data on My* 
presented in table 2 (for example, they are in the correct order of size) 
when allowance has been made, as in (4.14), for the fact that belly orbits 
do not give m* directly. In §6 we shall make comparison of m,* with 
an optical mass, m,*. 
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The main point here is that the observed electronic specific heat in 
the pure metals is quite consistent with our model. But when we study 
the effect of alloying we run into a thick fog of mystery. Almost all the 
experiments have been done on Cu, where, as Rayne (1956) pointed out, 
the reported values of y could be correlated with the estimated purity 
of the sample. For example, y for 99-999% Cu is 0-686—0-688, whereas 
samples of only 99° purity gave values as high as 0-75-0-8. ‘Thus there 
were either serious experimental errors or there is some effect which is 
extremely sensitive to small concentrations of impurity. 


Fig. 11 
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Variation of electronic specific heat with electron density for alloys of Cu with 
Zn and Ge (Rayne 1957 a, 1958). 


This has been to some extent confirmed by more systematic experiments. 
When Zn, or Ge, are added to Cu, we get the change of y shown in fig. 11 
(Rayne 1957a, 1958). Although there is a good deal of scatter, one can 
draw a reasonable curve through these points as if y depended only 
on the number of free electrons added with the divalent or tetravalent 
metal. On a free electron model with a spherical Fermi surface, a 
general increase of y with n is quite intelligible; the density of states 
increases as the Fermi surface swells until it touches the zone boundary 
in a special little cusp at about 1-32 electrons per atom. 

But another experiment spoils this argument. Guthrie (1959) has 
measured y for alloys of the type Cu,_,,Ni,,Zn,, which should retain n= 1. 
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For «=0-03, y increased from 0-686 to 0-844 (y/yp=1-68, and 
for «=0-08 y became 0-985 (y/yyp=1-98), which is much longer than 
any of the values attained in fig. 10. It seemed, also that the increase 
could be analysed into contributions as if from the Zn and Ni acting 
alone, so that it does not depend very much on the concentration of 
electrons. (On the other hand, in disordered Cu,;Au (Rayne 1957a), 
y= 0-66, in spite of a lattice expansion which would have increased the 
value for a free electron system to about 0-74, so that the effect of alloying 
may be here to decrease y). 

Our model, also, does not agree with the data of fig. 11. As we 
increase n, we expect to go down the curve, as we leave behind the little 
peak at the point of contact. As we see in fig. 8, the theoretical trend 
in y is always downwards, once we have passed the contact region. 

Cohen and Heine (1958) have suggested that the explanation of these 
discrepancies is to be found in a change of shape of the Fermi surface 
with alloying. In particular, they propose that the Fermi surface should 
get more spherical, even breaking away from contact, at least until so 
many electrons have been added that even a single sphere touches the 
zone boundary. But it is difficult to fit such a suggestion in our model. 
Let us suppose that the effect of alloying is simply to change the parameter 
u, the splitting at the band edge. Then looking at figs. 9 and 11, we see 
that we could only match the observed rise in y by making wu increase, 
from 0-2 say, to 0-4—surely an impossibly distorted surface. Sphericizing 
the surface—reducing w—does not increase y at all—it makes it fall even 
more steeply than the theoretical line drawn for u=constant in fig. 11. 
Thus, our numerical calculation does not at all support the interpretation 
put forward by Cohen and Heine; it suggests an exactly contrary effect. 

There is a further experimental point that may be of relevance. Biondi 
and Rayne (1959) have measured the shift of the main optical absorption 
edge in Cu (presumably due to transitions from the top of the d-band 
to the Fermi level) on the addition of Zn, and compared this with the 
theoretical shift, 


a» oi 
66 = aot ° . e ° e ° e (5.4) 


which can be calculated directly from the electronic specific heat data. 
There is a considerable discrepancy; the observed shift rises much more 
slowly than one calculates, and even when the curve becomes linear the 
slope is only about one third of the value given by (5.4). 

In the same experiments the absorption edge which is attributed to 
transitions across the energy gap into the upper band was less pronounced, 
and seemed to shift to longer wavelengths as Zn is added. This looks 
like the same sort of effect—the Fermi level rising towards the lowest 
state of the upper band. But our model tells us that once contact has 
taken place the energy gap for ‘vertical transitions’ in k-space remains 
constant, so that we have here some support for the Cohen and Heine 
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hypothesis that the band gap decreases on alloying. However the 
change is not big enough for contact to be lost (consult fig. 9 to see how 
u would vary with n in the series of surfaces containing more and more 
electrons but just keeping contact with the zone boundary) and would 
certainly not help with the problem of understanding the Hume—Rothery 
rules. The optical data and their theoretical interpretation are still too 
uncertain for this argument to be accepted without question. 

There are a number of points at which one might attack this problem. 
Is it possible that m* is radically changed—increased—with alloying ? 
There is no other evidence on this point. Again the alloying may cause 
a local distortion of the Fermi surface near the zone boundary increasing 
AN (&) there without having much effect on the round belly of the surface. 
One can show (Mr. W. G. Chambers, work in progress) that there could 
be different changes of energy at different points on the Fermi surface, 
due to scattering by impurities but the effect seems to be small. Then 
Jones (1957) has calculated an effect, that would look like an electronic 
specific heat, due to the change in the zero-point energy of the lattice 
waves as the electron distribution varies with temperature. If we thought 
that we were just about to make contact with the zone boundary, this 
would be important. Yet again, we know from the work of Pines (1955), 
that the electron—electron interactions and the plasma oscillations can 
seriously affect the density of states. It may be that the effect of alloying 
is to quench this term or to alter it in the way observed. Another 
possibility lies in Friedel’s discussion (1954) of the creation of virtual 
bound states in the neighbourhood of impurities. Finally, we may find 
that there is a flaw in the argument which derives the linear specific heat 
at low temperatures from the density of states in a model of independent 
electrons. 

This problem also merits further experimental study. For example, 
it would be interesting to check whether Ag and Au show the same 
anomalies and also to look at a wider range of alloys of Cu. 

A phenomenon which is related to the density of states is the Knight 
shift. Cohen and Heine (1958) have discussed this for alloys of the 
noble metals, but it is difficult to draw any firm conclusions from the 
measurement of a quantity which depends on a variety of complex factors. 
We should only like to make one comment. According to our model, the 
wave function in the state k is given by (3.1): 


y Py = Pye + Meg Peep? . . . . . . (5.5) 

where, as in (3.3), 
My g/% = (Fy — Oi)! Cag See te ernie) 
Thus, at the zone boundary itself, where dy-g=?_, and where 
E,-Ey= — |[%_,|, the state y, will be s-like or p-like according as 
U_,is —veor +ve. A naive interpretation of the argument of Phillips and 
Kleinman (1959) would suggest that the residual quasi-potential should be 
essentially attractive at large distances, so that we might expect Y%,,, to be 
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—ve, and these states always to be s-like—hence contributing to the 
Knight shift. However, a recent band-structure calculation by Segal 
(private communication) shows the p-state lowest in Cu, and some rather 
circumstantial evidence from the behaviour of the relaxation time (Ziman, 
to be published) makes it likely that this is the case also for Ag and Au. 
At this point the deficiencies of the quasi-potential model become apparent : 
we need proper band-structure calculations to argue in detail about the 
wave functions. 


§ 6. OpricaL Mass 


As pointed out by Cohen (1958), the optical properties of a metal 
depend in a special way on the shape of the Fermi surface, etc. Under 
appropriate conditions one finds that the electrons should behave as if 
they were free particles of mass 


ma = 12nhinN | | vd PD: sn GT) 


where there are N atoms per unit volume and the integral means an 
average of the electron velocity over the Fermi surface. 

It is not difficult to compute this integral for our 8-cone model. Thus, 
from (4.12), (5.1) and (5.2), 


* Za 
[ras= <= | v2 dz 
h ms 


= [Be 2fe)+tfeyla, ©. . (6.2) 


so that a/ fue ie 
(12) e/a (=e, 

rs 2 eee 

ma=minr/ sa 2ele 2) + 2u*In {r — 2) + +/{u? + (1 —25)?} 


= {tan <1) —tan-y="2)} |" ye GR) 


which can easily be calculated from values of z, and z, which have already 
been tabulated. In fig. 12 we plot m,/m* against n, for different values 
of u. We note that there is no apparent singularity in these curves at 
the point where the Fermi surface touches the zone boundary; the zero 
in v? cancels the infinity in the integrand of (5.1). In fig. 18, we also 
plot the optical mass as a function of wu, for a strictly monovalent metal. 

The experimental data of Schulz (1957) provide the entries in table 5: 
the optical mass for Cu is much larger than it is for Ag and Au. 
But we consult table 3 for appropriate estimates of w, and then read 
off, on fig. 13, theoretical values of m,/m*. Thus we have an estimate 
(m,*) of m*, the overall scale parameter of the energy surfaces. These 
estimates agree quite well with the values of m,* obtained by trying to 
match the specific heats; considering the inaccuracy of the experimental 
values of m,, this is satisfactory. 
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Optical mass, as a function of electron density, calculated for various valus of w. 


Fig. 13 
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Optical mass, as a function of w, for 1 electron per atom. 
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It is interesting to make this comparison in a slightly different way by 
calculating the ratio m,/m,=y/m,y, which ought to be independent of 
m*, but which should depend on the shape of the Fermi surface, etc. It 
turns out that this ratio does not depend so much on the parameter w 
as on the actual area of the contact region. One can plot (fig. 14), 


Table 5 


m, (obs) 


M,/m* (theor) 


hence m,* 
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Theoretical relation between optical mass and neck radius compared with 
experimental values. 
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m,|m, against r/p, the radius of the necks. Since we have values of 
r/p directly from the de Haas-van Alphen experiment, we can put an 
experimental point for each metal against this curve. The result is not 
very good, although the theoretical curve lies nearly within the quoted 
limits of experimental error. But this calculation certainly confirms the 
qualitative argument of Cohen for the conditions under which m,/m, 
should be greater or less than unity. 


§ 7. Hatt Errect 


Although a great deal of information about the Fermi surface can be 
obtained by the study of the effect of a strong magnetic field on the 
transport properties, this information is quite blurred and smudged at 
high temperatures and in alloys by the scattering of the electrons by 
lattice waves and impurities, and only the ‘low field’ phenomena can be 
observed except in very pure specimens at very low temperatures. Of 
these, the Hall effect is the strongest (being linear in #) and is relatively 
easy to measure, being a single scalar parameter in a cubic metal and 
hence well-defined for a polycrystalline specimen. 

For the theory we again refer to E & P (§12.5), where we find that 


e? ds 
R=o,/057; 09> male Sihaus eee a .(7.1)(7.2) 
and 
; e E x, Or dr \\ d® 
oy = ral {rit(vetMyy 5 =O Mar Fy — nat Oy 5 -».5-)} oR ‘ 
(7.3) 


In other words, oy is the electrical conductivity in zero field, governed by 
a relaxation time 7, which may vary over the Fermi surface. The Hall 
coefficient, R, obviously depends mainly on the components of the 
inverse mass tensor of the electrons on the Fermi surface, 
epgh “BOE 
) ye 12 0k,,0k,, 

A particular set of rectangular coordinate axes, with the field in the 
z-direction, is implied in (7.3), but because of the cubic symmetry the 
answer should be irrespective of their orientation. This symmetry 
allows us to drop all the terms containing derivatives of 7,, which vanish 
on the average over the Fermi surface. Then for a single cone we 
should choose the axis of the cone as z axis, whereupon M,,. | disappears ; 
the form (4.3) of &, is such that the cone axis is automatically a principal 
axis of the inverse mass tensor at all points on the surface region belonging 
to that cone. Thus, the value of o, is obtainable rather simply as the 
value of (7.3) for a single cone averaged for all possible orientations 
relative to the magnetic field, i.e. 


(7.4) 


= 2 2 = } . d 
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(7.5) 
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Unfortunately, we do not know 7,. But, if it were a constant, it would 
vanish from (7.1). Let us make this assumption for the moment. Then 
(7.5) can be evaluated by the same techniques as were used in (6,1)-(6.3). 


Fig. 15 
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Inverse Hall coefficient as a function of electron density, calculated 
for various values of w. 


We note from (4.3) that 
Meg t=Myytolfm*, Mi t=f"eyim* . . . (7.8) 
and by a little algebraic manipulation we obtain 


ee 2m in M21" (1—2,) 
= EE Sp ak AEE ea cia ARM EN LEAS es Fd Be Ys 
i nNec al m 34/3 te | \ (7.7) 


which can be computed quite simply for different values of u and for 
different energy surfaces. 
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We note at once that 

1 

ier is ee 
when m,=m* and u=0. This is, of course, the standard result for 
free carriers, the sign of R being determined by the sign of e. We see 
that R does not depend on m*, but on the shape and volume of the 
Fermi surface. It is conventional and convenient to re-interpret (7.7) 
in the spirit of (7.8), and to write 


(7.9) 


O O05 | (5 r/p *2 
Inverse Hall coefficient calculated as a function of neck radius. 


as if the Hall coefficient gave the effective number of free electrons per 
atom. This is the quantity which is plotted in fig. 15, against n, the 
actual number of electron states per atom available inside the Herel 
surface, for various values of w. 

There are several features of interest in this figure. The Hall coefficient 
is obviously rather sensitive to the shape of the Fermi surface; n* can 
be as large as 2, even for a surface with one electron per atom, for values 
of w in the range suggested in table 3. But we may also note that, on 
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each curve, n* is about 1-3 when the surface just touches the zone 
boundary. This suggests, as indicated in fig. 16, that n* depends mainly 
on the size of the necks, and not on the shape of the broad belly of the 
surface. We can plot a single curve which represents quite accurately a 
simple relation between n* and r/p, independently of the total volume 
inside the surface. 

What makes n* increase thus—or, more directly, what causes |R| to 
decrease as the Fermi surface expands? In semiconductor physics R 
is the balance of a contribution (conventionally negative) from electrons, 
and a contribution of opposite sign (i.e. positive) from ‘holes’. Thus, 
we should say that “the number of holes is increasing as the Fermi 
surface makes contact and pushes outwards over the zone boundaries”’. 
This sort of statement is often made, but it is slightly misleading in the 
present case. To achieve a true ‘hole’ carrier, we need to go right into 
a corner of a zone, where all the components of the effective mass tensor 
are negative. Only then can we construct wave-packets corresponding 
to the absence of an electron from a given state—wave-packets which 
behave like particles with positive charge and mass. The situation on 
the necks is more complicated. The components M,,-' and M,,~1 in 
(7.6) are certainly positive—i.e. ‘electron’-like—but the value of 
M,, 1 is negative. Thus, from (4.4), 


2 
f'(2=1- rarer - 7) asz—>1, een GlalO) 


and is negative for u<1 over the greater part of the range of z on the 
Fermi surface. (It is interesting to compare 1/f’(z), which is a sort of 
‘component of effective mass in the z-direction’, with the cyclotron-mass 
component (4.14). As we see in fig. 8, they have quite different 
properties.) That is to say, the effective mass tensor is hyperbolic, and 
a wave packet will obey rather strange dynamical laws, unlike those of a 
simple classical particle. It is true that we can distinguish between 
‘electron-like’ and ‘hole-like’ orbits in a strong magnetic field according 
to whether the orbits enclose full or empty regions in k-space, but a 
given electron state can participate in either type, according to the 
orientation of the magnetic field. Returning to (7.5), we see that the 
sign of the contributions to R will depend on the orientation of the 
electron velocity to the magnetic field. The thicker the necks, the 
greater can be the contribution from negative values of M,,-*, and hence 
the smaller the value of |R]. But that is not quite the same thing as 
saying that the number of ‘holes’ has increased. 

Now for the observed values of the Hall coefficient. Table 6 shows a 
comparison between theory and experiment. As shown in fig. 17, there 
is a good deal of uncertainty about the behaviour of & at low temperatures 
—probably because this parameter also is sensitive to impurities—but 
the room temperature data are fairly well defined (Onnes and Beckman 
1912, Wortman 1933, Chambers 1956, Frank 1957, Blue 1959). From 
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the observed value of R we calculate n* as in (7.9). For comparison we 
refer directly to fig. 16, using values of r/p from table 2 to read off 
appropriate theoretical estimates of n*. 

It is at once evident that there is a large discrepancy between theory 
and experiment; for all three metals, the observed value of n* is too 
small by something like 20%. The theoretical value is independent of 
m* go that we have not this parameter to play with, and the observations 
are quite accurate enough to rule out experimental error as the source 
of the discrepancy. 


Table 7 


Cu Ag Au 


Rops (em3/coulomb x 10~°) —4-9--— 5-2 —8-0+-—8:5 —7-0+—7-2 


Ry=(Nec)— (free electrons) — 7:45 — 10-65 — 10-60 


n* ops 1-52-+1-43 1-331-25 1-52-+1-48 


N*oale 2-05 1-58 1-87 


200 400 600 *K 800 


Variation of inverse Hall coefficient with temperature. References (F) Frank 
be (B) Berlincourt 1958; (L) Love 1959; (I) Fukuroi and Ikeda 
e Oe . 


One line of thought (Ziman 1959b) is to suppose that the thermal 
scattering of the electrons has to some extent sphericized the Fermi 
surface, so that at room temperatures the necks in Cu and Au are reduced 
to r/p~0-1, whilst in Ag contact has been lost with the zone boundary. 
This is what we should need if we wanted to obtain the observed values 
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of n* from fig. 16. But other evidence is against it. Thus, the measure- 
ments of optical mass ($6) made at room temperature are in good agree - 
ment with the specific heat data and the shape of Fermi surface found 
at low temperature by the ‘topological’ experiments. The change of 
Hall coefficient with temperature (fig. 17) is also inconsistent with this 
hypothesis, which would require n* to rise to the high calculated values 
of table 6 as we go down to low temperatures; we observe instead a low 


peak, and an irregular fall to smaller values in the residual resistance 
range. 


Fig. 18 


impurity |O atoms%O:lO free |-l electrons 2 per atoml-3 


Inverse Hall coefficient in alloys. Data taken from Andrewartha and Evans 
(1941), Blue (1959), Dorfman and Zhukova (1939), Flanagan and 
Averbach (1956), Frank (1955), Schindler and Pugh (1953); Wortman 
(1933). 


These irregularities in n* at low temperatures suggest that the Hall 
coefficient is sensitive to impurity type and concentration. When we 
actually study the effect of alloying, we find a bewildering confusion of 
phenomena. Figure 18, taken from various miscellaneous sources, 
shows how the room temperature Hall coefficient can be changed by 
10 or 20% by the addition of one or two atoms per cent of a polyvalent 
metal such as Sb or Ge. Moreover, this change is often in the wrong 
direction. Our theory would predict that n* should rise along one of 


D 
P.M.S. 
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the curves u=constant in fig. 15, as the electron concentration increases ; 
the typical behaviour is a steep drop, followed by a steady rise which is 
more or less parallel to the line for ‘free electrons’, but displaced from it 
by varying amounts for various types of impurity. On the side of the 
transition metals it is true, we find n* falling as it might due to electrons 
being swallowed by the d-shells, but the agreement with theory is only 
qualitative. 

Of course we may try to explain these phenomena also as the effects 
of ‘spericizing’ by the impurities. But the arguments of Cohen and 
Heine would suggest that each curve in fig. 17 should drop steadily down 
to the free electron line, and follow this up to the Hume-Rothery phase 
boundary at about 1-3 electrons per atom. There is no allowance in 
their theory for the necks in Cug)Si,9, say, being thicker than in pure 
Cu—which is what seems likely from fig. 16 and fig. 18. Moreover, the 
specific heat data of fig. 11 indicate a close similarity between CuZn 
and CuGe, when plotted in terms of electrons per atom; the effects of 
these two metals on the Hall coefficient are quite different. 

A more likely explanation is that put forward by Coles (1956) and by 
Cooper and Raimes (1959): the relaxation time 7, is by no means 
isotropic, and cannot simply be cancelled from (7.1), (7.2) and (7.5). 
In these expressions it occurs as a weight function—to the first power 
in (7.2), but squared in (7.5). Since our Fermi surface is highly aniso- 
tropic, there is a good reason to expect the electron scattering probability 
to be very different on the necks, say, from what it is on the bellies. 
Suppose, for example, that 7 is smaller near the zone boundary then it 
is on the main, nearly spherical part of the surface. The contributions 
of the neck regions in (7.5) will be reduced—and these are the parts 
where M,,~* has its large negative values which make |R| smaller than it 
would be for free electrons. Thus this anisotropy in 7 would have the 
effect of increasing the theoretical value of the Hall coefficient bringing 
it nearer to the simple free-electron value. That is, n* will be decreased 
towards unity. In crude terms, the ‘mobility’ of the ‘holes’ is less 
than that of the ‘electrons’, and they contribute proportionately less to 
the Hall effect. 

This is our interpretation of the discrepancy between theoretical and 
experimental values of n* in table 6; 7, for phonon scattering of 
electrons at room temperature must have just this sort of anisotropy, 
and our estimate of n* is to be reduced accordingly. Now, following 
Cooper and Raimes, we can imagine that the functional form of T,» 18 
different for impurity scattering from its form for phonon scattering. 
When impurities are added we first get a rapid change of n* as the residual 
resistance of the impurity becomes dominant over the ‘ideal’ resistance 
due to the phonons (one can check that the minima in fig. 18 came at 
about this point). Thereafter there is a steady rise in n* along a line 
characteristic of this form of 7,, due simply to the swelling of the Fermi 
surface as more and more electrons are added to the alloy. 
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Unfortunately it is difficult, as yet, to make a quantitative theory of 
this effect. The calculation of Cooper and Raimes is not applicable in 
detail, since it refers to a Fermi surface which is simply-connected and 
deformed only moderately from a sphere. But we can make a number 
of comments on the data. Evidently the degree of anisotropy of 7, 
is different for impurity scattering by different metals. It looks as if 


Fig. 19 


|.O I] 2 Ps gee! 


Variation of inverse Hall coefficient with electron density in CuZn alloys 
(Frank 1955). 


the anisotropy for scattering by Zn in Cu is actually less than for phonon 
scattering. This is confirmed by the curves of fig. 19 showing how n* 
is systematically less at 600°c than at 0°c, for the CuZn alloys ; at higher 
temperatures the phonon scattering is relatively more important, and 
tends to keep n* down. Generally speaking there is a tendency for the 
anisotropy to be greater (n* smaller) as the valence difference between 
solute and solvent increases; thus CuZn and AgCd lie above CuGe 
and AgSb. The observations on alloys of the noble metals with one 
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another are also most interesting. According to Onnes and Beckman 
(1912) in the AuAg system (fig. 20) where there are no disturbing changes 
of volume, n* rises at low temperatures to a high peak—above even the 
theoretical value already calculated for isotropic 7,. This peak is 
somewhat muted at room temperature, but is still quite remarkable, 


Fig. 20 


Au 20 40 60 80 “bo Ag 


Variation of inverse Hall coefficient Ye une AuAg system (Onnes and Beckman 
912). 


considering that the number of electrons is constant. We can only 
assume that the functional form of 7, for scattering by Ag in Au is quite 
different—is, say, larger on the necks than on the belly—from its shape 
for phonon scattering. However, the data of Flanagan and Averbach 
(1956) do not show such a large effect in AuAg, and in CuAu there appears 
a slight minimum in the plot of against concentration. Here is a point 
requiring careful experimental confirmation. 


Ordinary Transport Properties of the Noble Metals 35 


There is independent evidence for the anisotropy of 7 over the Fermi 
surface. Rayne (1959) has measured the infra-red absorptivity of Cu 
as a function of temperature, and in order to explain his results requires 
that 


CCl) Sel 04a tee fe age A Ns pol Tell) 


in which we mean that crys [and | {vd SS Se UB ©) 


and <1/r) is a similar average of 1/7, on the Fermi surface. Since this 
quantity is identically unity for isotropic relaxation time, we must 
conclude that 7, varies considerably over the surface—by a ratio of 5 
to 1 or more. (7.11) refers to electron-phonon scattering. Similar 
experiments on CuZn alloys (Biondi and Rayne 1959) give 


a) UO ihe eda oii a) sd eh) 


which suggests that the anisotropy of scattering by Zn is smaller than 
by phonons, in qualitative agreement with our deduction from the Hall 
effect. Although the theoretical deduction of (7.11) from the optical 
data is not quite complete, it would obviously be of great interest to do 
the same sort of experiment on alloys like AuAg, or like CuGe, where 
we expect the anisotropy of 7 to be entirely different. 


§ 8. MAGNETORESISTANCE 


The theory of the galvanomagnetic effects of second order is much 
more complicated, even for low magnetic fields where the method of Jones 
and Zener can be applied. Referring to E & P (§ 12.5) we find the follow- 
ing formulae: 

In a polycrystalline specimen there are only two independent magneto- 
resistance coefficients, the longitudinal and transverse coefficients in 
the relation 


Ap 
Po 


for the relative change of resistance in a wire at an angle 6 to the field H. 
These in turn are defined by 


EBSA ERIN IENG EU a RTS Sol) aes ete) 


B= p+ es Br=pi+ d+ ps? . . . . (8.2) 


in terms of the three quadratic-field coefficients for a cubic crystal, i.e. 
the coefficients in the formula 
fate : 
Po Po 


(where j,, h, are the direction cosines of the electric current and the 
magnetic field referred to the cube axes of the crystal). 
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Now the coefficients in (8.3) can be calculated from the following 
relations : 


2 
py? =— PoiF1102 + pol np |e 
po” = — po{O1111 — P1122 — P1212 — F201} a, 
2 
3 = — po{r212 + 1221 — PoTH' > 


where o,, is given by (7.5), and the other symbols are the components ofa 
tensor of the fourth rank: 
: et 7] 7] dS 
Capyd— — Epyv€uda 73,248 [reese {rege (rev} ars: 
(The symbols ¢,,,, ete. are, of course, the standard totally anti-symmetric 
tensors of the third rank arising out of the vector products between 
magnetic field and electron velocity in the equation for the Lorentz 
force.) 

These formulae are quite intractable unless we know the functional 
form of 7,. Suppose we make the simplifying assumptions of §7—that 
7 is constant and that the energy surfaces are of our model form. For 
the components of (8.5) that we need in (8.4) we find 


(8.5) 


(8.6) 


4111 = F212 = 9, 
e* > dS 
Cu 578 CoS = +M,,—1)+ ... cyclic perms} Ph (8.7) 


— 2e4 
1221 3 Bh 


: 1 
[0M 7M + ... cyclic perms} iis Primes t-8-V 


These expressions are relatively simple because, as previously remarked, 
the inverse mass tensor has the same principal axes over the whole of 
a cone, so that all off-diagonal components of M-+ can be made to 
vanish. 

Since all the terms in (8.7) and (8.8) can be expressed as algebraic 
functions of z, the integrations can all be done analytically. But this 
calculation would be much less realistic than for the Hall effect, because 
of the anisotropy of 7,. There is a very simple test. Add the relations 
(8.4): 

Pi™ + pa + pax™ = = pe* iiss «sn oe 
which by (8.6) should be zero. The experimental evidence (Olsen and 
Rodriguez 1957) is that this sum is not zero for any of the noble metals; 
but is comparable in magnitude with each of its separate terms. Thus, 
our model is not in agreement with experiment on an important general 
point. 

It may be that the shape of Fermi surface that we are studying is not 
near enough to the actual shape. A sufficient condition for (8.6) is that 


O. =fi(kz) +fo(ky) +f3(k,) a . $ : . (8.10) 


over the whole cap of a cone. This is only approximately true; as we 
go outwards from L towards K, say, in fig. 1, we shall notice the influence 
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of the contact region of the next hexagonal face, and this will tip the 
axes of the inverse mass tensor. Going towards U, on the other hand, we 
enter the realm of the more distant square face, which may not be so 
important. Thus, the surface loses the local rotational symmetry about 
each (111) axis, and there will be contributions to o,,,, from off-diagonal 
elements of M-! relative to this axis. 

A more likely explanation rests upon the anisotropy of 7,. Not only 
does this occur as a weight factor in the galvanomagnetic tensor (8.5) ; 
its derivatives in k-space should appear explicitly in the integrals. 
Looking at 0o,,,, for example, we find in the integrand terms like 


0? OT, \2 
rant nese + (5) | Cae eas al) 


which certainly need not vanish, even in cubic symmetry. If the 
anisotropy of 7 is really as large as has been suggested in §7 (e.g. (7.11) 
and (7.13)), then we must expect important contributions from such 
terms. The calculations of Olsen and Rodriguez (1957) are quite 
unreliable quantitatively in these circumstances, except in confirming 
that the Fermi surface touches the zone boundaries. 

This interpretation can be correlated with the observation of de Launay 
et al. (1959) that Kohler’s Rule is not exactly obeyed in Cu. When 
Ap/p) is plotted against H/p,, the points corresponding to different 
temperatures or specimens containing different impurities do not fall 
precisely on a single curve. The fundamental condition for the validity 
of Kohler’s Rule is that the functional form of 7, should be the same 
under all conditions, even though its magnitude may vary with 
temperature etc. We have already seen that this is unlikely from the 
Hall data; it seems that in Cu the anisotropy of 7, (e.g. the value of 
{r)<1/7)) is least for scattering by Zn, is larger for phonon scattering 
at room temperature and greatest for scattering by polyvalent 
impurities such as Ge. We might expect this variation to be reflected 
in deviations from Kohler’s rule. In particular, 6, the longitudinal 
magnetoresistance in a polycrystalline specimen, seems from (8.2) and 
(8.9) to be rather sensitive to the value at o,,,,, and hence to the form of 
7, It is interesting that B, is a bit smaller at room temperature than it 
igs at 78°K, as if the anisotropy of 7 had decreased and also that B, is 
not at all constant from specimen to specimen at 4:2°K, where the 
residual resistance probably arises from various atomic species as 
impurities. 


§ 9. ELEcTRICAL CONDUCTIVITY 


The theory of the basic transport coefficients is complicated, and it 
is difficult to say much quantitatively about the result of assumed 
deviations from a simple free-electron system. Thus, we need to know 
more about the electron-phonon interaction, more about the electron 
wave functions, more about the lattice spectrum, and more about how 
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to solve the Boltzman equation before we can assign definite effects to 
the anisotropy of the Fermi surface, or can calculate in detail the anisotropy 
of the relaxation time. The following remarks are not, therefore tied 
down quantitatively to our 8-cone model, but are concentrated on some 
points of difference between the observed conductivities and some of the 
conventional results for free-electron models. 

The magnitude of the ‘ideal’ electrical resistivity in the noble metals 
is not anomalous. When allowance has been made for the velocity of 
sound, the atomic weight, density, etc., of the solid, the reduced resistivity 
R (E & P, §9.7) is very similar in each case to the values noted for the 
alkali metals, and there is no systematic trend. The value of ® does 
not seem very sensitive to the shape of the Fermi surface, but there is 
no complete calculation to verify this theoretically. 

The temperature variation of p,; is enshrined in the parameter Op 
(E & P, $9.6), which turns out to be very close to the ordinary Debye 
temperature ©, derived from the lattice specific heat (see table 7). 
This in itself has some significance; it shows that the electrons are not 
being scattered entirely by longitudinal modes, whose Debye temperature, 
©, is much larger than Op, but that a great part of the electrical resistance, 
especially at low temperatures, is associated with scattering of the 
electrons by the shear waves of the lattice. 

We can look at this from two related points of view. At low tem- 
peratures we are talking about long-wave phonons, whose interaction 
with electrons can be represented by a ‘deformation potential’ (KE & P, 
§5.6). For reasons of symmetry, this interaction involves only the 
longitudinal modes if the Fermi surface is a sphere, but a multiply- 
connected surface should be much more sensitive to shear strains in the 
lattice, and there should be strong interactions with the transverse phonon 
modes, especially in the neck regions. 

The alternative view-point is based on an analytical theory of the 
electron—phonon interaction (J. G. Collins, to be published). In the case 
of free electrons the matrix element for the transition from state k to 
state k’ is given by the Bardeen method (E & P, § 5.7); 


I(k,k’)=e,°(k—k')E(|k—k']) ww. (9.1) 


where @(K) is approximately 3é, at K=0 but falls away rapidly as 
we go to large values of K. The phonon polarization vector is e,, so 
that if q=k’—k we find a factor e,*q which selects only the longitudinal 
lattice modes. To get interaction with the shear modes, we must study 
the electron—phonon U-processes, where 


o,2(K +k) aeon(q eee. ee oe) 


That is, when k’—k is too large to lie inside the Debye sphere we add a 
suitable reciprocal lattice vector, g, to q. Evidently (9.2) need not 
vanish automatically for transverse modes in which e, is normal to q. 
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Now if we study the same matrix element for transitions between 
wave functions like (3.1), we get a more complicated formula; one can 
easily verify by substituting that this takes the form 


eg" | he" VA, dr = Sota eM(kKtgk’ +g) . . (9.3) 
&8 


where ¥% (k, k’) is much the same as in (9.1) and the coefficients « are 
the solutions of the secular equations (3.3). It is at once obvious that 
(9.3) contains terms with polarization factors like ene (Kea a= Kes g); 
which will give contributions from the transverse modes even when 
k’— k=q. Thus, there is no need to invoke U-processes in order to get 
interactions with the shear waves of the lattice. 

Here it is worth making a short digression on the subject of nomen- 
clature. In the free electron model we can label states by their actual 
wave vectors, so that we know the absolute values of k’ and k, and can 
verify that their difference, as in (9.2), is larger than any possible q, 
and hence requires an Umklapp process for scattering. Much of the 
discussion of the theory of electrical and thermal resistance is concerned 
with the different roles played by N-processes and by U-processes—as 
we have just noted, we make statements such as ‘transverse modes take 
part only in U-processes’, etc. But in the repeated zone scheme this 
distinction is only conventional. We can choose to say that the wave 
vector of our final state is 

Me Ke Cer ert Gaus act |e a (Oe) 


(as in fig. 21) because electron wave-vectors are arbitrary up to the 
addition of a reciprocal lattice vector. We then have 


CK Kk wh sea eens...) (9:0) 


and the transition has the formal properties of an N-process. Of course 
the polarization factor is still (9.2) (because this depends on the actual 
wave vector of the free-electron state) but the classification scheme has 
begun to fail. 

When we have a multiply connected Fermi surface, as in fig. 22, the 
distinction between U and non-U becomes entirely artificial. Suppose 
that k” lies on a neck, just across the zone boundary from k. Are we to 
distinguish the transition to such a state from the transition to k, say, 
which lies on the same neck but just inside the zone boundary. In 
the repeated zone scheme the energy surfaces are continuous along the 
necks, and the zone boundary is only a conventional plane for counting 
unit cells in reciprocal space. It is more natural to describe the transition 
by the shortest possible vector, from k to k”, than to add on a reciprocal 
lattice vector, go to k’, and then have to subtract the same reciprocal 
lattice vector from k’—k in order to find the phonon wave vector q. 
This seems to be the best convention to use in (9.3) say, trusting to the 
values of the coefficients «, etc., to give the right contributions from 
the various OPW’s. One can show, for example, that the matrix 


40 J. M. Ziman on the 


element is a continuous function as the final state vector crosses the zone 
boundary, so that there is a smooth change-over from * V ~processes ’ 
to ‘U-processes’, without any abrupt alteration in the physical 
phenomena. 

In this formulation, can we understand the: anisotropy of 7, The 
problem is very complex, and can only be solved numerically (Ziman, to 
be published). There are a number of different competing factors, such 
as the following : 


Fig. 21 
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N-processes and U-processes on a multiply-connected Fermi surface. 


In the first place, the coefficients a, ,, ete., vary considerably over 
the Fermi surface; according to (5.5) and (5.6) we have something like 
a single OPW over the belly regions, but this is modified into a mixed 
state, with two OPW’s of nearly equal amplitude as we go to the zone 
boundary. If %,,, is negative all the terms in (9.3) will be positive, so 
that there will be an enhancement of the scattering matrix element in 
the neck regions. 
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Then we must look at the density-of-states factor in the scattering 
probability. If the bulk of the scattering is through small distances on 
the Fermi surface (and there is a tendency for this, especially at low 
temperatures) a high local density of states, as at a neck, will give a high 
transition rate, and a small value of 7,. 

A third factor lies in the definition of 7, itself. This is by no means 
a simple question, for we are really studying the solution of a complex 
integral equation and there is no guarantee that there is a solution even 
of the form (cf. E & P, §§7.7, 9.5) 


eareV EER Ie) ny (ete uy) ae (9.6) 


in the electric field E. But if this is so one can argue, as a rough 
approximation, that 


1 ; 
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where 9," dk’ is the probability of a transition from state k into a 
state k’ in the range dk’ and (v,, v,,) is the angle between the electron 
velocities in those two states. Now, if 9," has a strong forward 
lobe, the value of the integral will be larger on the necks than on the 
belly because the radius of curvature is smaller on the necks, so that 
cos (V,, V,-) deviates more quickly from unity as we take k’ away from 
k, and the weight function has less effect in suppressing the contribution 
of ‘local’ scattering to the integral. This argument is equivalent to 
the statement we make in the free electron model where a U-process, 
of quite small phonon wave vector, can reverse the direction of the electron, 
and hence make a large contribution to the electrical resistance. (cf. 
fig. 21). 

Can we also understand why the anisotropy of 7, should vary from one 
scattering mechanism to another? This is a much more uncertain 
problem, but there are arguments above which make 7, depend on the 
shape of the differential scattering cross section of the impurity or phonon 
displacement, and we may expect this to vary a lot when we compare 
the scattering by the screened Coulomb field of a charged impurity 
(e.g. Blatt 1957) with the effect of a homovalent impurity such as Ag in 
Au. This effect, also, has been confirmed by a detailed calculation 
(Ziman, to be published), and the behaviour of 7, found to be consistent 
with the Hall effect in alloys. 

How is it then, that the conventional theories of the residual resistance 
produced by chemical impurities work so well? As shown in E & P 
(§9.2) there is good agreement between the experiments and a theory 
of free electrons scattered by screened charged impurities. Until we 
have exact solutions for the scattering matrix elements of the proper wave 
functions, this is an open question, but it may be simply that the relaxa- 
tion time of the electrons in the neck is so short that they play little part 
in the conductivity and the belly states carry the current and are scattered 
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nearly as if free. The attempt to calculate by the variational method 
(Ziman 1959b) may be seriously misleading, for it may greatly over- 
emphasise the contribution of the strong scattering in the contact regions. 

Whilst discussing electrical resistance we naturally think of the 
‘resistance minimum’ phenomenon, and the corresponding anomaly in 
the thermoelectric properties. It now seems (Gold et al. 1960) that this 
effect is always associated with the presence of free Fe, or other magnetic 
impurities, and is not related directly to the shape of the Fermi surface. 
An interpretation in terms of breaking away of contact with the zone 
boundary (Ziman 1959b), is not supported by the other evidence, and 
must now be abandoned. 

Linde (1958) has observed changes in the apparent slopes of the 
resistance/temperature curves in Cu and Au when alloyed with certain 
other metals, and has noted a correlation with the pressure coefficient 
of the additional resistivity. This is most likely to be due to changes 
in the Fermi level (on thermal expansion or elastic compression of the 
lattice), which sweeps through the supposed broad ‘resonances’ in the 
scattering cross-section of transition metal impurities (Friedel 1956). 


§ 10. THERMAL ConDUCTIVITY 


At high temperatures this is given (presumably; has it ever been 
thoroughly checked?) by the Wiedemann—Franz law, and introduces 
no new results. The low temperature value is of some interest because 
there are certain theoretical comparisons that can be made (Klemens 
1956). Thus, one can calculate 


D;=1/%,=64T?W ,(0)/O2W,(T) . . . . (10.1) 
D,=1/Y9,=0-12902p,(T)/Ly>T3W,(T) . . . (10.2) 
(for notation see EK & P, §9.10) which are tabulated in table 7. Despite 
some uncertainty in the proper choice of © (we simply follow Klemens) 
there is obviously a significant difference between the observed values of 


these parameters for the noble metals and the theoretical Bloch-type 
free-electron model which would have D,;=D,=1. 


Table 7 


Cu Ag Au 


Op (°K) 310-335 | 210-225 | 165-186 


Or (°K) | 320-333 | 200-223 | 170-200 
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We notice at once that the magnitudes of D, and D, are correlated 
with the diameters of the necks on the Fermi surfaces. Is this significant ? 
The conventional argument is that these ratios are sensitive to the 
degree to which U-processes dominate the electron scattering, especially 
at low temperatures; there is an important difference between electrical 
and thermal conduction (and between thermal conduction at high and 
at low temperatures) in that the former is particularly sensitive to 
scattering through large angles, as in a typical U-process right across 
the electron sphere, whereas heat conduction at low temperatures is 
affected chiefly by the loss of energy of an electron when it collides with 
a phonon, irrespective of the angle of scatter. For a nearly spherical 
surface one can discuss these ratios in terms of the shape of the electron— 
phonon interaction (E & P, $9.10), and it is easy to show that D, and 
D, need neither be equal to unity nor to one another (Cohen and Heine 
discuss the role of different averages of sound velocity in D, and D,, 
but this is not the whole story by any means). For our more complex 
surfaces it is difficult to draw any quantitative conclusions from the values 
in table 7. We can note that we have the equivalent of U-processes in 
the scattering of electrons across regions of high curvature on the Fermi 
surface, and these should, according to (9.7), have a larger effect on 
p, than on W,, but it does not seem possible as yet to compute D, or 
D, directly and to see how they should depend in detail on the shape of 
the Fermi surface. 

Indeed there is a difficulty about (10.1) and (10.2). It is assumed 
that 7T?/W,(7') is a constant at low temperatures—that the ‘ideal’ 
thermal conductivity obeys the theoretical Bloch—Wilson law xoc7'-*. In 
the case of Cu there is strong experimental evidence against this; Powell 
et al. (1959) favour a function like 7-8, which would make D, and D, 
into variables increasing as the temperature falls. But really this 
should not surprise us. The 7’ formula is based upon very special 
assumptions about the form of the electron—phonon interaction. It can 
be shown (Ziman 1954) that this formula must be modified when we 
allow for U-processes, even on a nearly spherical Fermi surface. For a 
surface with contact, and with a complicated interaction with phonons, 
we might get considerable deviations from a simple quadratic law. By 
the same tokens, the 7° law for the electrical conductivity is not 
sacrosanct; deviations from it have been reported (Powell et al. 1959) 
although the full experimental position is not yet clear. 

Another phenomenon showing how far we are from the conventional 
Bloch theory is the deviation from Matthiessen’s rule, for both electrical 
and thermal resistivity in cold-worked Cu (Powell et al. 1959). But 
this is just what we would expect from our hypothesis that the anisotropy 
of the relaxation time of the electrons is quite large, and varies from one 
type of scattering to another. It can easily be shown (EK & P78 7-10) 
that there must be positive deviations from Matthiessen’s law—the 
resistivity of a system containing several species of scatter must be larger 
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than the sum of the resistances due to each species separately—if 7, 
is different for the different types of scattering. It is dangerous to combine 
data from different sources, but it is worth noting that Berlincourt (1958) 
found the Hall coefficient of cold-worked Cu to depend considerably on 
the temperature, as if the anisotropy of 7, for scattering by dislocations 
etc. were very different from the anisotropy due to phonon scattering at 
200°K. It would be interesting to make a systematic search for deviations 
from Matthiessen’s rule in specimens showing unusual values of the Hall 
coefficient, and see whether these effects are correlated. 


§ 11. THERMOELECTRIC PoWER: THE ‘ DiFFUSION ’ TERM 


By tradition, the thermoelectric properties of metals and alloys offer 
important clues to their electron structure. For example, it is well 
established theoretically that a simple free-electron system should have 
a negative thermopower—just what we observe at room temperature in 
all the alkali metals except Li. So, when the thermopower is positive, 
we feel that the Fermi surface must be complicated, with a prepon- 
derance of ‘hole’ states etc. The fact that @ is positive for the pure 
noble metals then seems strong qualitative evidence for the sort of 
multiply-connected surfaces which we are here studying. 

But when we look quantitatively at this evidence, it is not so impressive. 
On the contrary, the theory of the thermopower of the noble metals is 
as uncertain as many of the experimental results which it sets out to 
explain! 

First the observations: fig. 23 shows the temperature dependence of 
Q, the absolute thermopower of pure Cu, Ag and Au, at low tempera- 
tures (Pearson, private communication) and up to the melting point 
(Cusack and Kendall 1958). These data are probably the best at present 
available but they have not the weight of being agreed by a number of 
independent authorities. The points to notice are (a) that Q is positive 
except at very low temperatures, (b) that there is a hump in Q at low 
temperatures, (c) that for Cu it varies linearly with absolute temperature 
above about 200°K, (d) that the line for Ag does not go through the origin 
and (Gold, private communication) tends to bend over above 1000°K, (e) 
that Au does not show a definite linear portion above room temperature, 
but rises to a broad maximum and then falls slightly. 


The standard theory gives a very simple formula for the electronic 
contribution; thus (K & P, §9.11) 


TR dlno(é) 
Set eee an EE 
Qe= T= | ws ‘bigs eke ioc eel Lees 


where o(&) is ‘the electrical conductivity if the Fermi level were at &’. 
In this formula e is the conventional charge of an electron, that is, a 
negative quantity. For a simple free-electron system, we write 


o(S)ocr(P)uy(P)H(P) . . . . . . (11.2) 
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showing how the conductivity depends on the Fermi energy through the 
relaxation time, the Fermi velocity, and the area of the Fermi surface. 
It follows very easily that 


Q mRkRT (3 o{In 7(&)} 

ee (ee Ne | 

e003 6€25\2 O{In &} 
7 R kT 


where ¢ is a dimensionless constant, which, for phonon scattering at 
high temperatures should lie between +2-5 and +3. 


Fig. 23 


200 400 600 IO °K ICCD 1200 


The thermoelectric power of the noble metals (Gold e¢ al. 1960, Pearson, 
private communication, Cusack and Kendall 1958). 


This theory tells us that Q, should be directly proportional to 7’. 
If we look at the curves of fig. 23 below 300°K, and assume that the 
hump is due to phonon drag and has disappeared at about 200°K, we 
find for € the values —1-7 for Cu, —1-1 for Ag and —1-5 for Au. But 
these linear trends are not continued at high temperatures except for Cu. 
For Ag the slope above 600°K corresponds to €~ — 2:5, whilst for Au the 
value of é is quite uncertain. Thus, the experimental parameters are not 


well defined. 
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Nevertheless, we want to explain these negative values of €. Is the 
funny shape of the Fermi surface a sufficient cause? Let us look at the 
rate of change with energy of the analogue of v,© 1.€. 


fi=6,| 2{Infoasl [ac | asa eas 


This has the value 3 in (11.3). From (6.2), using (4.6) and (4.7) we find 
a simple formula for ¢,, for our 8-cone model 


é =( 1 )( 2 i. 2(2,—21) +f (2) — (S42 + Lf (ex) PY Ser +f (ea)} 
Lien Ne AN eae PP ‘ 
eo Tak | (2e—2f(2) + Lf (Pde 


: (11.5) 
In fig. 24 this quantity (with m*=1) is plotted against the number of 
electrons per atom, for several values of wu. We see that it falls steeply 
from about 1-5 as the Fermi surface expands, and there is a change of 
slope at the contact with the zone boundary. Yet €,, remains obstinately 
positive until the surface holds 1.1 electrons or more, and certainly does 
not reach negative values of the order of —1. In other words, if the 
main part of the diffusion thermopower were simply due to changes in 
the mean Fermi velocity, and Fermi area, we should expect Q, to be 
between —0-5 and —Il1,pv/°c at room temperature, instead of the 
observed values which lie between + 1-5 and +2,v/°c. 

This argument is not, of course, complete, for it ignores the variation 
of + with energy. To allow for this, Bailyn (private communication) 
has suggested that 7(&) simply mirrors /vdS, so that €=2,. This 
obviously does not fit the observations; the argument for it depends on 
a crude use of the variational method, which is not really sufficiently 
accurate for such a delicate calculation. Jones (1955) has pointed out 
that for phonon scattering one might expect 7(&)océ&?/4(&), putting 
the emphasis on the probability of scattering being proportional to the 
density of states into which the electron can go. But if we look at fig. 9 
we see that this does not help to make Q change sign. It is true that 

Olnr(é) _ OlnW(é) 


Srlny tert ee . . . . . (11.6) 


can become negative as we approach the zone boundary—indeed ./(&) 
has infinite slope at that point—but beyond the point of contact the 
density of states falls again, and the contribution of (11.6) to € will be 
firmly positive. 

We may try to blame our model of the energy surfaces—to argue, for 
example, that the proximity of other zone boundaries may already have 
taken us beyond the point where 

Or, 106, /06, 

0é& hoe] ok ~ * 
has become negative over a large part of the Fermi surface (ef. fig. 8). 
But remember that the Hall coefficient also contains terms of this sort 
(M,,-* in (7.5)) and we cannot afford to improve our value of Q at the 


(11.7) 
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expense of R (which should already be smaller, theoretically, than we 
actually observe). In vulgar parlance, we have already fond the Hall 
effect showing more ‘electrons’ than ‘holes’ in the noble metals: now 
the thermopower suggests that ‘holes’ predominate over ‘electrons’. 


Fig. 24 
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Calculated variation of €,, with electron density, in the neighbourhood of 
contact. 


The only hope lies, again, in a more complete theory of 7(@). We 
must write 


a(6)o | ands ete teta Shi ttee a LHS) 


and calculate the derivative of the whole function. There will be effects 
due to 7 as a weight function in expressions like (11.4), but, more par- 
ticularly, there will be the explicit variation of 7, with energy. Going 
back to §9, we may expect large contributions from the variation of 
the coefficients «,, etce., with energy (cf. (5.6)). Our hypothesis that 
7, is smaller on the necks than on the belly of the surface is consistent 


P.M.S. E 
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with the need for 7() to decrease rapidly as & increases, i.e. as we approach 
the zone boundary in general. It is possible that we could find a link 
between a calculation of this effect and the magnetoresistance coefficients, 
as in (8.11). But all this waits on a good computation of the ‘ideal’ 
resistance of the noble metals sensitive to the details of the electron— 
phonon interaction. 

There is a considerable literature on the thermoelectric properties of 
alloys of Cu, and a few results for alloys of the other noble metals (e.g. 
Domenicali and Otter 1954, Domenicali 1958, Friedel 1956). It is 
generally observed that any element added to Cu reduces the thermo- 
power sharply—presumably because dIn7(¢)/dé is smaller for scattering 
by impurity atoms than it is for phonon scattering. But most of the 
analysis that has been made of these data has been based upon the simple 
free-electron model with scattering by screened, charges, etc. Now that 
we know the full complexity of the Fermi surface, it is extremely difficult 
to accept the results of this analysis, even though some of its qualitative 
features (e.g. ‘resonance’ scattering by d-shells of transition metals) may 
be retained. All that one can hope is that the anisotropy of 7, is so 
great that there is not much contribution from the neck regions, so that 
the thermoelectric and other effects come mainly from the belly, where 
the electrons are more or less free anyway. 


§ 12. THERMOELECTRIC POWER: PHONON DRAG 


It is now recognized that there may be an important contribution to 
the thermoelectric power from the heat carried by the lattice waves 
which are emitted as a current of electrons is drawn through the solid. 
This effect has been shown to be especially important at low temperatures ; 
in the alkali metals (MacDonald ef al. 1958, Ziman 1959a) the sign and 
magnitude of the ‘lattice thermopower’ may be correlated with the shape 
of the Fermi surface. We must therefore consider this phenomenon in 
the case of the noble metals. 

The formal theory (K & P, §§9.13, 10.9) is complicated; in the end 
we can use one or other of the following formulae : 


Listes Ch sit 
5 (ave) (pe): op pee new kath se 


V,= Ky e in 3Nk oP 
o,.'R cha ( Pu i; hed te ele (12.2) 


Here Cis the lattice specific heat, x, is the lattice conductivity, and o 
is the ‘ideal’ electrical conductivity. These formulae are in principle, 
equivalent, but the first is more convenient at low temperatures (say 
7’ <0-2@) and the second at higher temperatures. 

For the moment let us ignore the symobols P,,, Peli. end 
concentrate on the other factors, which fix the general scale and tempera- 
ture dependence of Q,. Remembering that k/e=—86 LV/°O, we see 
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that @, can be quite large. At low temperatures, where (12.1) is 
appropriate, the 7° specific heat law holds: 

sve = (5) op olan Be, Ga) 
Thus @,, could be of the order of 1 zv/°c at 7’ = 0/20, and should increase 
rapidly in magnitude above this temperature. But then the lattice 
conduction begins to be influenced by phonon-phonon interactions, etc. 


Lattice thermopower, estimated from fig. 23. 


(which makes it difficult to estimate P,,) and we go over to (12.2). As 
shown in E & P, §§ 8.2, 8.9, when we get above the Debye temperature, 


we can expect to write 


Salita ye Wea lee oo ORS) 
R\o,1' Poye VL) en 


where y is the Gruneisen parameter—say y~2. This again, if—P,,/Py; 
is not too small, means quite a large thermopower—perhaps 8,v/°c 
at T=@:; at lower temperatures it would rise exponentially as the 
U -processes are frozen out, until electron-phonon interaction dominates x,. 
E2 
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It is true that in Cu and Ag this would be reduced by isotope scattering 
of phonons, but a short calculation (KE & P, $8.6) suggests that 
this would be a negligible effect above room temperature. In fig. 23 we 
see well-developed humps in the thermopower at low temperatures which 
Blatt and Kropschot (1960) have shown to be due to phonon drag. If we 
suppose that @, is added on to a linear diffusion thermopower, we 
get the curves of fig. 25 (which are plotted in terms of 7/0, to make 
them comparable). ‘There is no doubt that the phenomenon is much 
the same in all three metals—Q, is positive, rising rapidly to a peak of 
about Lyv/°c near 7'=0-150 and a slower fall to zero at about 7'=0. 
Cu differs slightly by having a narrower peak shifted to about 0-20, 
and has, it seems, a small negative part at very low temperatures. 

Do our formulae (12.1)-(12.4) explain these phenomena? Having 
said nothing, so far, about (—P,,/P,,) and (- P,,/Py), let us treat 
these as adjustable parameters; we find that we may represent the 
general shape of the experimental curves by using (12.1) and (12.3), 
with (—P,,/P,,)~ —9:14 below the peak, and using (12.2) and (12.4), 
with (— P,,/P,,) ~ — 0-025 at higher temperatures. The major discrepancy 
is that the observed value of @, seems to fall to zero more rapidly than 
this model would allow at temperatures above 0-5@ (indeed it may 
become slightly negative, although this is difficult to check against 
the linear diffusion term). 

So we are forced to discuss these variational integrals P,,, P,,, P,,. 
In the theory of E & P these are defined as follows: 


Pu= pr | | | (0, — D, PP,’dk dq dk’, 


1 ; F 
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(strictly speaking P,, should contain further terms corresponding to 
phonon-phonon interactions, scattering of phonons by impurities, etc., 
but these are all taken care of by using (12.2) at high temperatures—i.e. 
in the range where scattering of electrons by phonons is no longer the sole 
agent governing «,). 

In E & P we further postulated that 


Dek tu; Oe q sul ces eee ee eee 


whence we find, for electron-phonon N -processes, that P,,= —P, tp=Pr- 
This would probably be the asymptotic situation at low temperatures 
for a spherical Fermi surface. But for our more complicated system this 
assumption is too crude. Fortunately (12.6) is not the essential condition 
for (12.1), (12.2) and (12.5) to yield an acceptable variational estimate 


Ordinary Transport Properties of the Noble Metals 51 


of @,. We can use any trial functions ®, and ®,, provided that they 


are normalized to make 
(fe)... aan 
z | er, Pi (Afi /96y) dk : 


4 | vaPa(Ong"/O2) dq 

ay e 
Naturally we should like to choose these trial functions so that they 
conform to the functional behaviour of the electron and phonon distribu- 
tions that actually satisfy the Boltzmann equation. We can do this by 
writing 
@, = Key -u; ®,=“%q-u “if lind ae Ns GW) 
h Ne UND, 
where 7,, T, are relaxation times for electrons in state k, phonons in 
state q, and 7,, 7, are average relaxation times for electrons and phonons. 
—i.e. defined so that the kinetic formulae 

Pean Nein Kp 4C,0,°7, 7... (12,9) 

give the actual ideal electrical conductivity and lattice heat conductivity. 
Of course we do not know 7, and r, directly, and can scarcely calculate 
them but at least we have some intuitive ideas about how they should 
behave. 

When we put (12.8) into (12.5), we get some complicated integrals for 
P,,, P,, and P,, which can be reduced a bit by the use of the methods 
of KE & P, §9.5, but which cannot easily be evaluated. P,, and P,, 
come out as positive definite functions, not necessarily equal. But the 
_ sign and strength of Q, will depend on —P,, about which we can make 
the following qualitative comments. 

At low temperatures, only long-wave phonons, with small values of 
q, are allowed. All processes in our formulation (cf. fig. 22) are V -processes, 
so that k’=k-+q, and we can approximate to the key factor in — P,,: 


e 
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Or : 
; (4 . st) 744). (12.10) 


The second term will tend to be small, because q will lie nearly in an 
energy surface, normal to v,. In other words, we see again the inverse 
mass tensor. If we suppose, even more crudely, that 7, is cancelled, 
in the integration, by the scattering probability factor, we are led to 
something like an average value of the components of M~ tangential to. 
the Fermi surface. To get some idea of how this might behave, we 
calculate 


ma |M, eM et Meret = /s [e-+ 45 (12.11) 
Zo &y4 


This is now our (wildly approximate) estimate of the limiting value of 


=mr,74(4°M~*+q)+ 


52 J. M. Ziman on the 


—P,,/P,,. Using the numerical data already amassed, one finds a 
number that is small, and positive, being about +0-1 for the sort of 
surface that we assume for Ag, a bit larger for Au, and perhaps + 0-2 
for Cu. If one looks closely at fig. 25, at the very low temperature end, 
one can say that the ‘observed’ values run from about —0-3 for Ag, 
through —0-07 for Au, to a small positive value for Cu. These are in 
the right order, but shifted to more negative values. Obviously we 
cannot expect numerical agreement, but we do see that there is no 
inconsistency between the thicker necks in Cu and a slightly negative 
value of Q, at the lowest temperatures. The simple statement that 
‘when the surfaces contact the zone boundary U-processes are possible 
down to the smallest phonon frequencies, and these make @, positive’ 
is only roughly correct. 

At higher temperatures, when larger values of ¢ are allowed, it becomes 
impossible to discuss these integrals in terms of the local differential 
geometry of the Fermi surface. Above the Debye temperature we have 
transitions between all parts of the surface, weighted by the electron— 
phonon matrix element, the density of states, etc. In the conventional 
theory, these transitions will contribute positively to —P,,/P,, (ie. 
negatively to the thermopower) if they are V-processes, but may subtract 
(i.e. give positive contributions to @,) if they are U-processes. In our 
model, whence U-processes have been banished by law, we need another 
criterion. From the fact that —P,, contains the factor 

G* (TM — THY) 

this comes down, more or less, to the statement that the contribution to 
Q,, will be negative or positive according to whether the chord q between 
the points k and k’ passes through occupied or unoccupied regions of 
the Fermi surface (fig. 26). In fact this criterion is similar to the rule 
for calculating the sign of the contribution of an ‘orbit’ to the Hall 
effect in a strong magnetic field; if the orbit encloses filled states in 
the zone, then it is an ‘electron’ orbit; if it encloses empty regions, 
then it is a ‘hole’ orbit. Transitions between points on an ‘electron’ 
orbit then give negative lattice thermopower; transitions across a 
“hole” orbit give a contribution corresponding to carriers of positive 
charge. This argument is not quite exact if the Fermi surface has a 
very complicated geometry, but it does avoid the language of 
“N-processes’ and ‘ U-processes’. 

Unfortunately, without a long and complicated calculation, it is 
impossible to verify the experimental conclusion that (Py; 1P 7) 
small and negative in the three noble metals—although perhaps it goes 
to zero around 7’= ©. No explanation is offered for the extremely similar 
values of @ ar the maximum in all three metals, although we may note 
that this peak value is rather less than the maxima in Rb and Cs, which 
go up to 3pv/"c. (MacDonald et al 1958). There is conflicting evidence 
(Powell et al. 1959, Pearson 1955, von Oijen 1957) on the effect of cold 
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work on the thermopower of Cu. It now seems (Pearson 1960) that these 
effects are due almost entirely to traces of Fe impurity, but it would be 
interesting to see whether, by changing the dislocation concentration and 
also altering the isotropic constitution of the Cu (c.f. Blatt and Kropschot 
1960), one could change @, through a change in the form of T, Finally, 
it must be agreed that there is little hope of explaining the Eel 
behaviour of the thermopower of Ag'and Au at high temperatures (fig. 23) 
by appeal to phonon drag, which by (12.4) will surely have disappeared 
above the Debye temperature. 


Fig. 26 


“electron” transition | 


“hole” transition 


§ 13. Some UNDOUBTED Facts anD UNCERTAIN THEORIES 


It is not easy to sum up this complicated subject. The noble metals 
have suffered so many experimental observations, which have been 
subjected to so many theoretical explanations, that one can always find 
some evidence to support any point of view. Many of the observations 
are inconsistent with one another, and with allthe theories. The abnormal 
sensitivity of many properties to small amounts of impurity leads to 
conflict of evidence, which only further deliberate experiments can 
resolve. It is not merely a question of fitting together a jig-saw puzzle ; 
every piece has several spurious versions which must be identified and 
discarded. 

But the overall conclusion is simple; there is no serious inconsistency 
between the observed transport properties and the assumption that the 
energy surfaces discovered by the ‘topological’ techniques remain more 
or less rigid when the temperature is raised or when the metal is alloyed 
with other metals. If there is any sphericizing effect of alloying, it does 
not seem enough to draw the Fermi surface out of contact with the 
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zone boundary. It is a pity that we cannot, therefore, maintain the 
well-known and elegant explanation of the Hume-Rothery rules—that 
the phase boundary in alloy systems comes at a definite electron-atom 
ratio corresponding to contact of a spherical Fermi surface with the face 
of the zone—for now we must believe that there is already substantial 
contact in all the noble metals. Perhaps another geometrical criterion 
will be found to explain this phenomenon, see Hume-Rothery and Roaf 
(1961). 
More detailed points that have been discussed are as follows: 


(i) The band gap at the zone boundary should lie in the range of 
5-10ev to be consistent with the shape of Fermi surface found by the 
topological techniques. 

(ii) The change in the electronic specific heat with alloying cannot be 
explained by a simple alteration of the shape of the Fermi surface. 

(iii) The observed optical mass is consistent with the shape of the 
Fermi surface and the electronic specific heat. 


(iv) The Hall effect cannot be fitted without assuming considerable. 
anisotropy of the relaxation time of the electrons. 


(v) The change in the Hall coefficient with alloying is probably due to. 
differences in the anisotropy of the relaxation time for different scattering 
mechanisms. 


(vi) The low-field magnetoresistance depends very much on the aniso- 
tropy of 7, and cannot be used for estimating the shape of the Fermi 
surface. 


(vii) Deviations from the simple Bloch—Wilson theory of electrical and 
thermal conductivity are qualitatively explicable, but cannot yet be 
interpreted quantitatively, although there is good experimental and 
theoretical evidence for fairly strong interaction between the electrons. 
and transversely polorized phonons. 


(viii) Substantial anisotropy of the relaxation time can be explained 
from the anisotropy of the Fermi surface and variation of the electron 
wave functions in the matrix element of the electron-phonon interaction. 

(ix) The positive thermoelectric power at high temperatures is not due 
to changes in Fermi velocity or Fermi area, but to a rapid decrease of 
7(é) aS we approach the zone boundary. 

(x) Phonon drag is the cause of the low-temperature hump in the 


thermopower, but this does not provide much quantitative information 
about the shape of the Fermi surface. 
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§ 1. THe Nearty INDEPENDENT QUASI-PARTICLE MopeEL 


Tue usual electron theory of metals developed since the 1930's has been 
built up on the independent-electron model, i.e. on the assumption that 
the electrons can be treated as if they did not interact with each other 
at all. However the coulomb repulsion between the electrons is a very 
strong force, and the purpose of this paper is to discuss the electron theory 
of metals taking the coulomb force into account from the beginning. 
We shall not discuss many-electron effects such as plasma oscillations, but 
limit ourselves to transport phenomena. For the purposes of this paper, 
the term ‘ transport phenomena’ will be taken to cover all genuine transport 
properties such as electrical and thermal conductivities, as well as other 
properties like the electronic specific heat and the magnetic spin suscepti- 
bility that depend only on the electrons at or near the Fermi level. The 
thesis we shall develop is that such properties can be understood in terms. 
of quasi-particle excitations at the Fermi level, which to a first approxi- 
mation are independent so that one can still use the language and frame- 
work of the independent-electron model: but the short-range screened 
coulomb interaction existing between the excitations leads to some 
important corrections and modifications of various types which we shall 
describe in detail. Indeed, this paper explores how far it is or is not. 
possible to cast the exact theory of the electron gas into the preconceived 
mould of an independent-electron model or at least a nearly-independent 
quasi-particle model; for the independent electron model has served 
solid-state physics so well for so long and will probably continue to be 
used for a long time in analysing and correlating the various electronic 
properties of metals. 


1.1. Introduction 


In the independent-electron approach, each electron is in a Bloch 
state #,(r) with energy ¢«(k), and all states with energy «(k) less than the 
Fermi energy ¢,, are occupied by electrons. (We shall use ¢ for the energy 
of single particles and reserve EF for the energy of the whole system of 
electrons.) The ground state of the electron gas is a single determinant, 
of all the #,(r) with ¢(k) <e,, corresponding to the assumption that each 
electron moves completely independently of all the other electrons. The 
Fermi surface defined by ¢(k)=e, is an immediate natural concept, 
dividing the occupied and unoccupied regions in k-space. The excited 
states of the system are still single determinants, but now with some of 
the Bloch functions #, with k=k,Jk,....... from below the Fermi surface 
omitted, and an equal number with k = k;, kj, .... from above the Fermi 
surface included. As is often done, we shall describe such a state by 
ignoring the unchanged part of the Fermi distribution and simply saying 
that we have some excited electrons with k=k,, k ‘eres and some holes. 
with k=k,,k,.... (We shall continue to distinguish electrons and holes 
by Roman and Greek letter subscripts.) Such states can be excited. 
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thermally and by electromagnetic fields, and it is on this independent- 
electron model that all the usual electron theory of metals found in the 
standard texts is based (Mott and Jones 1936, Seitz 1940, Wilson 1953). 
As is universally known, the independent-electron theory has served 
solid-state physics extraordinarily well, and it is still the basis of almost 
all discussions of transport phenomena; see for instance Ziman (1960). 
At first sight this is rather surprising, because the electrons in a metal 
cannot be moving independently of one another, and the coulomb repulsion 
between them is in fact a very strong force. For instance, in sodium the 
electrostatic energy of two electrons at the interatomic distance apart 
is 4ev while the Fermi band width is only 3ev. Thus the electrons 
accelerate and decelerate one another, and even in the ground state there 
is ample energy available to excite some electrons temporarily into high- 
momentum states. The wave function for the ground state of the system 
is very complicated, and if we expand it as a sum of determinants of 
Bloch functions, it will certainly contain some %, with energies well above 
the Fermi level. Thus we might expect the Fermi surface and all simple 
properties associated with it to be more or less completely washed out. 
While this is true in the simple sense of the independent-electron model, 
there are strong experimental and theoretical reasons for believing that 
a sharp Fermi surface in some more sophisticated sense does exist and 
that much of the independent-electron model can be salvaged provided 
it is interpreted in a modified way. As one piece of experimental evidence, 
consider for example the de Haas—van Alphen effect (Shoenberg 1957), 
which indicates the existence of very small pieces of Fermi surface enclosing 
in some metals less than 10~° electrons per atom, forming pockets perhaps 
at corners of the Brillouin zone. The radius of the pockets is about 2°, 
of the radius of the Fermi surface, and the maximum variation of «(k) 
within the pockets about 0-04ev. Also most de Haas—van Alphen oscilla- 
tions are unobservably weak due to thermal excitation at already a few 
degrees absolute where &7’~10ev. As Mott (1956) has emphasized, 
such results reveal some very fine detail in the electron distribution, and 
would be quite inexplicable if everything were smoothed out by lev or 
more and no sharp Fermi surface existed in some sense. Another argument 
is due to Skinner (1940) who observed sharp emission and absorption 
edges corresponding to the Fermi level in the soft x-ray spectra of metals. 
From the discussion of the last paragraph, one might expect a broadened 
edge due to the electrons temporarily in high-momentum states. However 
the point to note is that the soft x-ray spectrum corresponds to the energy 
levels of the whole system, the one-electron fluctuations inherent in the 
ground state being irrelevant. Consider a state of the system in which 
one electron has been given some extra energy « above the ground state. 
It will rapidly collide with other electrons, exciting them and sharing its 
energy with them. Since the extra energy any one electron can get varies 
between 0 and «, we would expect that the number of different configura- 
tions that the original state can decay into would be at least proportional 
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to «. Thus as we consider a state with excess energy tending to zero, we 
would expect the lifetime + of the excitation to tend to infinity because 
there are no states for it to decay into. This would imply that the states 
of the whole system with infinitesimal excitation energy have zero energy 
broadening Aexhi/7, leading to a sharp edge in the soft X-ray spectra in 
agreement with experiment. It also implies that the excitation can move a 
long distance through the lattice, and experimentally mean free paths of 
the order of 1mm have been reported for very pure specimens of tin and 
other metals at low temperatures. One can define a wave vector k for 
the excited state in the usual way: 


Wir, +R,.tatR,,----)=exp Gk, Ro yoe (05 tos nes) eee ree eee 


where “ is the complete wave function for the state and R, any lattice 
displacement. The existence of a sharp energy in the limit of infinite 
lifetime, coupled with a k vector, again suggests a sharp Fermi surface 
outside which k for the excitation has to lie. Finally we may consider a 
semiconductor like germanium where the mean electron—electron electro- 
static energy (several ev) is larger than the band gap (about 1 ev) and at. 
first sight quite able to excite electrons across it. Indeed the ground state 
of the system will contain high momentum fluctuations as before, but we 
know experimentally that these cannot carry a current and contribute to. 
transport phenomena, because with high purity material at low tempera- 
tures free carrier concentrations of 10-1! per atom can be achieved. 
However, if we put an extra electron into the crystal, it can move about, 
having strong electrostatic interactions with all the valence electrons, 
locally exciting them and temporarily being excited by them which will 
involve high-momentum components in the wave function. We can now 
picture the whole mess of electron plus the local disturbance it causes. 
moving together through the lattice with well-defined energy and wave 
vector and a long free path leading again to a one-electron-like band. 
structure ¢(k). 

This is about as far as the experimental evidence itself can take us— 
not at all conclusive but nevertheless suggestive of the existence of a 
Fermi surface, sharp in some sense at least to about 10-tev, with low- 
energy one-electron-like excitations having well-defined wave vector and 
energy and a relatively long lifetime. 

On the theoretical side, much sophisticated work has been done 
recently, and this is reviewed in §2. Although a complete and rigorous 
solution to the problem is not yet possible, nevertheless a fairly clear 
picture does emerge from the theory to date. This picture, which we shall 
now describe, is in complete accord with what has already been con- 
jectured from the experimental side and indeed puts it on a much more 
detailed footing. Indeed Landau (1956) realized early on that a picture 
of this sort must apply, and analysed some properties of the electron gas 
in terms of it. A very readable comment on the whole situation has been 
given by Herring (1960). 
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1.2. The Model Described 


Let the ground state of a gas of N electrons be denoted by Yow or |GN > 
in the Dirac notation or |G) for short, and its energy by ty ee arena Cle 
very complicated wave function, embodying all the correlations beween 
the movements of the electrons and their accelerations and decelerations. 
We can always define a total k vector for a many-electron state using (1.1), 
and we assume that V has been adjusted to the symmetry of the system 
so that the ground state is non-degenerate with total k=0, e.g. N has to 
be even. 

Now there is a minimum energy ¢,, the electron affinity of the system, 
required to add one extra electron to the system. If we add the electron 
with an energy « greater than ¢,, we can produce certain one-excited- 
electron-like states |k;) of the V+1 particle system, which Hugenholtz 
(1958) has shown can be expressed in the form 

i ees Cy eters etn pet se dC) 
where the creation operator, O,;* is a definite operator defined by 
eqn. (2.31) below. Since the electron interacts with all the other electrons, 
we cannot ascribe to it an individual Bioch function and wave vector, 
but it follows from general principles of group theory (see for example 
Heine 1960) that the lattice symmetry of the crystal and (1.1) can be 
used to define a total k-vector k; for the state as a whole. We shall in 
what follows take the symbol k; to include also a specification of the spin 
direction. The expectation value of the energy for the state (which in 
general is not an eigenstate, see below) we shall write as 


ieee UK) ger saden ebay, | oles yin, oe) (1.3) 
and we shall describe «(k;) as the energy of the electron-like excitation 
which we have put into the system. For large values of e«(k)—e, our 
whole description breaks down, but for energies less than an ev or so 
the wave vector k; always lies outside a precise surface defined by 

(ier ere martes ub eoreae so (14) 
This, the Fermi surface, has the same size as in the independent-electron 
model, i.e. it contains exactly one electron per atom or whatever the 
valency of the metal is. (For a proof of these and other statements, see 
the references reviewed in § 2.) 

The minimum energy required to remove an electron from |G, ) is 
also ey. This is one of the characteristics of a metal, for it is equal to the 
energy that has to be added in going back from |G, VN —1) to |G, and 
this is the same as in going from |G, NV) to |@, N+1) by continuity. Only 
in an insulator or semi-conductor with a band gap at the Fermi level 
is there a discontinuity at NV electrons when the valence band is completely 
full and the next electron has to go into the conduction band. Now by 
extracting an electron from below the Fermi level, we can create certain 
excited states |k,) of the N —1 particle system 


Fe ONC) tices Cia cre eee aoe (10) 
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where the annihilation operator O,, is the Hermitian vonjugate of the O,* 
defined by (2.31). The expectation value of the energy of the state is 


B=H,~e(k),2 2 ee ee 


and we say the system contains one hole-like excitation of energy 
—e(k,). The function ¢(k,) is continuous with ¢(k;) at the Fermi level 
€,, and k, has to lie inside the Fermi surface for the state (1.5) to exist. 
Now the states defined by (1.2) and (1.5) are not in fact eigenstates of 
the system, and they decay with a lifetime 7'(k) into more complicated 
states (see § 1.3.1 and § 2.3). However 7'(k) is proportional to [e(k) — «,]~ 
so that the lifetime tends to infinity at the Fermi level. In particular the 
excitations produced thermally and by electromagnetic fields have a 
sufficiently low energy, that their lifetime 7'(k) is much longer than their 
relaxation time due to collisions with lattice vibrations and impurities. 


Fig. 1 


Multiple excitation corresponding to a uniform current. |@—electron-like 
excitations ; OQ—hole-like excitations. 


Thus we shall describe the excited states of the N-electron gas in terms 
of multiple excitations. 


[Ks Keys «a's Kez, Weps es) Og Og SO Ope «| ae Be 


with equal numbers of electron-like excitations (wave vectors Ky Kane 
and hole-like excitations (k,,k,,...). It does not matter that these states 
are not eigenstates, although they actually tend to eigenstates as all k;, k, 
tend to the Fermi surface and the life time becomes infinite. For finite 
energy of excitation, the proper eigenstates are very complicated combi- 
nations of all states of type (1.7) with the same total energy and same 
total k= >k;—>k,, and they do not seem to correspond to anything of 
physical significance. We can describe all normal excitations of the 
system in terms of the states (1.7). For example, a uniform electric field 
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produces electron-like excitations on one side of the Fermi surface and 
hole-like excitations on the other (fig. 1), equivalent to the usual uniform 
displacement of the whole Fermi distribution in the independent -electron 
model. 

Thus an excited state of the system can be represented as a low density 
gas of quasi-particles, namely electron-like and hole-like excitations. 
Although the density of electrons is of order one per atom, the density 
of quasi-particles is much less, e.g. of order &T'/(e,—¢))~ 10-2 per atom 
at room temperature where ¢,,—€,) is the band width, corresponding to 
the number of excited electrons in the independent-electron model. 
From the way in which the quasi-particles have been defined from (1.2) 
and (1.5), we expect them to be closely related to electrons and holes. 
Indeed the operator O,* and the properties of the quasi-particles are 
derived using perturbation theory from the corresponding independent- 
electron states by mathematically ‘switching on’ the coulomb-repulsion. 
However a quasi-particle differs from an actual electron (or hole) in 
several respects. For instance, it is not a point particle because it includes 
all the local disturbance in the electron gas which screens the charge. 
The screening results in the interaction between two quasi-particles 
being modified from a pure coulomb repulsion to a short-range ‘screened 
coulomb’ force. For present purposes we shall take it to have the form 


FenU OX DE Terma. ke ae. £138) 
where the screening distance 1/( is of the order of the interatomic spacing : 
a discussion of its precise shape is given in § 2.4. Finally, if we form a wave- 
packet out of our excitations, the charge e* transported by the packet is 
not exactly that of one electron. The whole question of the velocity, 
charge and current carried by the excitations is discussed in § 3. 


1.3. The Consequences of Interactions 


So far we have reviewed the individual properties of the quasi-particles, 
treating them as independent. We must now take account of the screened - 
coulomb interaction between them, and also discuss their lifetime 7'(k) 
which is a related effect. 


1.3.1. Auger transitions and the lifetime T'(k) 


Consider a state with a single excitation with wave vector k, and energy 
€,=«(k,). As already mentioned this is not an eigenstate of the system 
because it can make Auger transitions to other states with the same total 
energy and total k, the only two quantities that are proper constants of 
the motion (Skinner 1904, Landsberg 1949, Hugenholtz 1957, Quinn and 
Ferrell 1958). In particular because of the coulomb force, the excitation 
can “collide with an electron at k, below the Fermi surface’, be scattered 
to k, loosing energy ¢,—«,, and ‘eject’ an ‘electron’ into the excited 
state k, leaving a hole at k, (fig. 2). A better description would perhaps 
be to say that the original excitation has interacted with the Fermi sea 


P.M.S. F 
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and decayed into another excitation at k, plus an electron-hole-like pair 
at k,, k,. We must have 


k, =k, +k,—k, and €,=€,+ €3— €4. eS NS a 
It is easy to do an approximate calculation. The lifetime 7’ (k) is given by 
y) 
12 eh) ee 
fhe 
Fig. 2 
h z 
+] 


a ae 


(a) 
e hs 
() 
e h, 
(b) 
(a) Initial and (6) final states for Auger transition. (c) Construction for density 


of states. 


where AN (E) is the density of states for the whole system and M the 
matrix element of the Hamiltonian connecting the two states. We approxi- 
mate to M using independent-electron wave functions with plane waves 
and the interaction (1.8) : 
exp(—itk,:r,) exp(—ik,-r,) e? 
a a a manera 1 armies TNE, 
exp(tk;-r,) exp(ik,: ro) 
cir — qu arias 
47re” 
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where q=k,—k,=k,—k;. To calculate (EH), we note that k, can lie 
anywhere in the shell between energy ¢, and ¢,, i.e. can range over a volume 
47k,*A in k-space where A is the distance of k, from the Fermi surface. 
Having fixed ¢, and k,, we next fix «, which can have any energy between 
€, and ¢,—(<,—€,), ie. the surface «=e, can cut a radius of the Fermi 
sphere anywhere along a section of length about $A in k-space if we assume 
on the average that half of the energy is given to the electron-hole pair. 
k,; must now lie on the energy surface ¢,=«€,+(€,—«,). In order to satisfy 
(1.9), we redraw the surface ¢(k) = «, displaced by wave vector — q (shown 
dotted in fig. 2). k, can then lie on the circle of intersection of this surface: 
with ¢(k)=eg, ie. along a line of length 27($k,.) say, which then fixes. 
k, at k;+q. In order to get the density of states per unit energy, we allow 
k, to range from energy «, to energy «,+dH, giving a strip of width 
dH |(d«/dk). Multiplying by the density Q(27)-* of points in k-space, we: 
obtain 
MN (LE) & Q?(277)-* (4k, A) EAT, (0€/0k)—. 


Putting A=(e—,)/(0«/0k), q~Q~ky=1-9/rs, 0e/Ok ~ 2€,,/ky, substituting: 
in (1.10) and dropping all numerical factors, we obtain finally 


h e2\2 ] € 2 : 
aia el a) reer 


where rg is the radius of the sphere containing one electron. As stated in 
§ 1.2, n/T varies as (e«— €,)" giving a very long life for excitations near the: 
Fermi level: e.g. for «-—e,=871=3x10ev at room temperature, we 
obtain a mean free path of ~10-3cm which is much longer than that for 
phonon collision and therefore negligible. However (1.12) also shows. 
that for energies of more than a few ev, the lifetime broadening is sufficiently 
large for the concept of a band structure with single identifiable excitations. 
to be no longer very meaningful. 


1.3.2. Inter-particle collisions 


The Auger process discussed above represents the decay of a single: 
excitation by itself. At a finite temperature when there are many quasi- 
particles, there are also collisions between them. We can calculate the 
mean free path approximately as follows. For static obstacles of concen- 
tration c, a quasi-particle has to travel through 1/c atomic cells before: 
striking an impurity and we have L~rs/c. In our case the concentration. 
of scatterers, namely other quasi-particles, is of order &7'/e,. Also the 
momentum transferred to the struck quasi-particle has to be such that; 
the latter can conserve energy, i.e. the available k-space into which it: 
can be scattered is only a fraction &7'/e, of the whole space available., 
which reduces the scattering by this factor : 


L li! (BT ey). 


This formula is very similar to (1.12). Indeed the two processes are: 
very similar, collision with an unexcited electron just below the Fermi 


F2 
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level in the one case and with an excited electron or hole in the other, and 
they are usually treated together in making proper calculations. Such 
calculations are described by Abrahams (1954) and by Ziman (1960), who 
also reviews the relevant experimental data. The results briefly are these. 
In a free electron gas, electron-electron scattering cannot affect the 
electrical resistivity because the total k and hence the total current are 
conserved in each collision. However, with a lattice potential and 
complicated Fermi surface, it contributes a term AZ? to the electrical 
resistance both by Umklapp scattering and by scattering from a part 
of the Fermi surface with high velocity to one with low velocity. Electron— 
electron scattering also contributes to the thermal resistance even for a 
free electron gas, because a state in which ‘hot’ quasi-particles are flowing 
against a counter stream of ‘cold’ ones can clearly be spoilt by collisions 
between them. Numerically the effects seem to be just a little too small 
to be clearly detectable except for a 7 term in the electrical resistance of 
transition metals (Ziman 1960, Baber 1937, Mendelssohn 1956). 

In conclusion, the interaction effects of §1.3.1 and this section just 
give another scattering mechanism, whose contributions to the thermal 
and electrical resistances can be calculated by the usual independent- 
electron transport theory and anyway are in almost all cases small. Thus 
they do not give any new effect of an essentially non-independent particle 
nature, and we shall not consider them further. 


1.3.3. Breakdown of the exclusion principle 


Although the operators O,* of (1.2), (1.5) are too complicated for making 
explicit calculations with as yet, we shall assume that the exclusion principle 
applies exactly at the Fermi level and the operators O,* satisfy the usual 
Fermi commutation rules exactly there, since at the Fermi level single 
excitations represent eigenstates of the system. Thus one can only have 
one excitation of a particular k, and a state O,*O,*|@) is identically zero. 
This however may not apply at higher energies where the excitations are 
not eigenstates, and our whole mode of labelling and counting states 
begins to break down. For example, if we start with one electron gas at 
a low temperature and another at a high temperature and accelerate the 
electrons with electric fields in the same way, the resulting currents may 
not be exactly equal. Any deviations would presumbaly be proportional 
to (e—e,)?, giving a contribution in 7? to the electrical resistance. We 
shall assume that the effects are at most comparable with those of electron 
electron scattering discussed in the previous section, and not consider 
them further. 


1.3.4. Hnergy corrections to transport phenomena 


The interaction between quasi-particles gives corrections to the energy 
which in principle are important for all transport properties. If we 
consider as state of the system with two quasi-particle excitations k, and 
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k,, then the energy # of the system is not just the sum of their individual 
energies but includes a term 7(k,,k,) due to their interaction : 


Ei(ky, ky) = + €(k,) + e(K5) + (ky, k,). tee en (Le L8) 
(Here as elsewhere the minus signs have to be used for holes; compare 
(1.3) and (1.6).) Following (1.12) we write it in the form 


dees 
Ki, Ka) —A(k,,. ke) — Seer sige: ate (1914) 
n(k,, ky) (ky, a) 7 F ; (1.14) 


: 
where the last factor is a few ev and the first a numerical factor whose 
value we might expect to depend roughly as follows on spin and on whether 
we have electron-like or hole-like excitations: 


both holes or both electrons parallel spin Aw~ —0:2 
one hole and one electron parallel spin Aw~ +0-2 
both holes or both electrons antiparallel spin A ~ —0-05 
one hole and one electron antiparallel spin = A~ +0-05 


These are just rough values for ordinary metallic densities calculated from 
the theory of Bohm and Pines (Pines 1955, eqns. (6.4), (6.23)): in detail 4 
of course varies with k,, k, through ¢(k) and the electron density. 

In an ordinary transport process, multiple excitations with n quasi- 
particles are excited where n is of order V but a very small fraction of it : 
n/N ~&T/v,, or vp/vyp, Where vy is the drift velocity e&/mz and v, the Fermi 
velocity. The total energy is 

H=H,+#, 


=[De(k,)— Dek, )]+42 > 0(k, kk). . . . 2. (115) 
The two-particle interaction energy H, includes the double summation 
over all pairs of quasi-particles, and is therefore of order n?/N. As regards 
the one-particle energy H,, if the excitations were all by a finite amount 
Ae about the Fermi level, then #, would be of order » and the interaction 
energy negligible by comparison with it. But the multiple excitations one 
always has in transport phenomena are really distortions of the Fermi 
surface, the change in energy being quadratic of order n?/NV, so that the 
interaction energy is comparable with the one-particle contribution H,. 
As regards order of magnitude, #, can of course be larger than #, in metals 
like nickel or cobalt with small de/0k, where it leads to ferromagnetism. 
In ordinary metals it is perhaps 40% of #, if a change of spin and hence 
the exchange energy are involved, as evidenced by the corrections to the 
Pauli spin susceptibility (Pines 1955). When spin changes are not 
involved, corrections of the order of 10°, might be expected. This is 
probably a maximum estimate actually, because with a simple band 
structure one has a sum of an equal number of positive and negative terms 
on account of the equal numbers of electrons and holes. We take up the 
topic of the size of the corrections again in § 5, and shall see that the above 
figures based on the Bohm and Pines theory are seriously overestimated. 
The actual effective electron—electron interaction is bot weaker and more 
smoothly varying with k than the interaction of Bohm and Pines, resulting 
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in substantially smaller corrections. Incidentally the fact that the 
density of excitations is always so small, allows us to neglect three-body 
and higher interaction terms in (1.15). 

Although it is clear that the interaction energy affects quantities like the 
spin susceptibility which depend essentially on an energy balance, it is 
not immediately obvious whether it comes into genuine transport pheno- 
mena like the d.c. conductivity. Here the normal mode of calculation 
determines the amount of momentum given to the electrons by the electric 
field and hence the current carried, without explicit mention of the energy. 
But fundamentally electrical resistance is concerned with the energy fed 
into the electron gas and its rate of dissipation into heat, and we have to 
develop a slightly unorthodox method of calculating resistivity to bring 
this out. 

Every transport phenomenon is characterized by the type of multiple 
excitation it produces: in the case of d.c. conduction there are electron- 
like excitations on one side of the Fermi surface and hole-like ones on 
the other (fig. 1). The whole is equivalent to a distortion of the Fermi 
surface measured by some parameter A; in the present case A=dk/k, 
where 5k is the side-ways shift of the surface. The total energy of the 
excitations is 

Bea, © 535.0 See 


where the constant « contains significant contributions from the inter- 
action energy since like #, it varies as A?, as already discussed. The 
electric current carried by the system is discussed in §3: in any case it is 
proportional to A: 

J=BA. Lot hae eee aes 


The electrons are scattered by impurities and phonons, so that without an 
electric field to maintain it the current would decay with relaxation time 7: 


dJ /dt = — J |r. ne te eS ei cy 
Thus from (1.17) we can write the energy as 
Hf = (af B*)I*®), io er CR ee 
and its rate of dissipation as 
dE dJ 
ie ee) Cay dates 
dt (toe dt 
= 2a J? 
wee 


But the rate of dissipation of energy is RJ? where R is the resistivity, 
and we have 
Bix= 20, Bree nee (1.20) 
which shows explicitly how the interaction energy contributes to R 
through the value of a. 
In §4 we shall apply such arguments to some particular transport 
processes. 
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1.3.5. Current interactions 


It is in principle possible that the interaction between quasi-particles 
leads to corrections, analogous to 7 in (1.13), to the current carried by 
them. We can regard the n—1 other excitations as constituting a modi- 
fication of the background medium in which a particular quasi-particle 
moves, producing a change of order v,n/N in its velocity v,. Now the 
total current carried is of order nv,, i.e. linear in n (unlike the one-particle 
energy H, of eqn. (1.15)). In comparison with this, the total interaction 
corrections to the current are of order v,n?/N and can therefore be 
neglected. 


§ 2. THEORETICAL FOUNDATIONS 


In this section we review the most important contributions to the 
theory of the electron gas made in the last five years. We emphasize 
those works which, directly or indirectly, point to the existence of a sharp 
Fermi surface. In particular we shall be interested in the existence and 
properties of the low-energy excited states, defined in terms of quasi- 
particles. In § 2.1 we discuss the various methods of calculation and the 
difficulties that arise in evaluating the resulting expressions. In §2.2 
we review the theory of the ground state in the free electron approximation 
while § 2.3 deals with the excited states and the quasi-particles. In § 2.4 
the screened interaction between quasi-particles is discussed. In § 2.5 
we summarize briefly what has been done to include the periodic potential 
of the lattice, and in 2.6 the work on the derived properties of the electron 
gas such as specific heat etc. 


2.1. The Theory of the Electron Gas 


We first define our notation. The many-body hamiltonian for our V- 
particle system is 


H=H +A, . . . e e e 5 (2.1) 
where 
KH =k > ka, *a, + » » U(G)a,* ay 46 eer Pe ooo) 
k G+0 


gives the kinetic energy and the lattice potential energy of the electrons, 
and 


*y * 9 
HA ,=4 » Vx 23444 Ons Ugg . . ° . ° (2.3) 
1234 


expresses the electron—electron interaction energy. 

In this formulation a,* and a, are respectively the creation and anni- 
hilation operators of electrons of momentum k and spin o; G is a reciprocal 
lattice vector and 


U(6) =" [U(rhexp (iG nae Sram tera o*d.) 
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is a Fourier coefficient of the lattice potential. The matrix element 
Vy934 i8 defined by 


1 ‘ 
V1934 = Gk — ka)d(or, 54)8(09, o3)0(K, + Kg, k; +k), eee 
and - 
4ir/k? k40 26 
v(k) = { 0 Net ON (2.6) 


is the Fourier transform of the coulomb interaction. We have chosen for 
simplicity atomic units, i.e. 
AK=e=m=l1, 
and Q is the volume of the system (assumed to be indefinitely large). 
The eigenvectors of (2.1) must satisfy the conditiont 


Ya,*a,[PS= NID © a ee 
k 


expressing the constancy of the number V of electrons where the density 
of electrons py= N/Q remains finite as Qo. The interelectronic spacing 
is defined by the parameter 


3 1/3 
6s= ( ) . . . . . . . . (2. 8) 
tipo 


The solution of the problem consists now of finding the eigenvectors 
and eigenvalues of (2.1) which satisfy (2.7). In most cases the free electron 
approximation has been used; i.e. the potential due to the crystal lattice 
is replaced by a uniform background of positive charge, and in (2.2) we put 

U(G)=0: Bisa 
This drastic approximation greatly simplifies the calculation, and it is 
not definitely clear yet whether the results obtained for the free electron 
gas can be directly generalized for the lattice case. It is in general assumed 
that qualitatively this generalization is possible, although some important 
corrections might be perhaps necessary (see below). i 
To solve the Schroedinger equation 


AVE yeh|Y > 
use of perturbation theory has almost invariably been made, the inter- 
action #, being considered the perturbation term. A comprehensive 


discussion on the perturbation series by Brueckner (1958) shows that due 
to the large volume of the system (QQ-> 00) 


(a) the Brillouin—Wigner perturbation method gives an essentially non- 
convergent series and cannot be used to solve the problem of a many- 
body system ; 

(6) the Rayleigh-Schroedinger perturbation theory, although it seems 
to give a divergent series, can be used to perform the calculation because 
it can be shown explicitly that the divergent terms, through an exact 


(Aa Ree RE, Sit a a le de ae 
+ We shall always use angular bras and kets for perturbed states, and round 
ones for unperturbed. 
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cancellation of their contributions, give no net contribution to the final 
result of the energy value. All the perturbation methods used to solve 
the many-body problem are modifications, in one way or another, of the 
Rayleigh—Schroedinger series. 

Let us firstly study the unperturbed hamiltonian #%,. Its eigenstates 
which satisfy (2.7) are vectors |’) such that N different k-states are 
occupied by electrons, the others being empty; these vectors |‘t’) are 
equivalent to the Slater determinants of the real space representation. 
The unperturbed ground state |G) consists of N electrons occupying the 
states 


k= |k|<k,, 
where £,, is defined by 
Or 1/3 
ky = (377? py)? = (Z) rs = 1:91973+ et lenne (2eLo) 
and its unperturbed energy is 
Wie ee eI at (2a 1) 


in atomic units. The occupation function f(k) is a constant equal to 1 for 
k<k, and zero for k>k,; it has thus a discontinuity Af=1 at k=ky,. 
The low-energy excited states of the unperturbed system can be defined 
as the ground state, plus n electrons in given states k,(k,> k,) and n empty 
states or holes k,(k,<k,). These excited states will be denoted by 
[k,, k;,...5 Ky, Ky...) Ae BS rhage 8 (22 
and their energy is 
WM KK hbe, Ken Kapaa) Went Wo 4) he 8 5 (2,13) 


In this way it is possible to redefine the unperturbed system in terms 
of unperturbed quasi-particles. In this new picture the ground state is 
considered the vacuum, and each excited state (2.12) consists of n electron- 
like quasi-particles of momentum k,, spin o; and energy $k,?, together with 
n hole-like quasi-particles of momentum (—k,), spin (—o,) and energy 
( = xk”). 

To study the perturbation we use the time-dependent adiabatic formula- 
tion. This method has been employed by Goldstone (1957) and Hubbard 
(1957, 1948 a). The other methods used by other investigators are slightly 
different. We do not attempt to describe them but shall give their results 
expressed systematically in terms of the present adiabatic formulation 
as far as possible. 

We must mention separately the works by Kohn and Luttinger (1960) 
and Luttinger and Ward (1960), where using a thermodynamic approach, 
they have tested the validity of the adiabatic approximation. In fact 
they have shown that Goldstone’s method of calculation is valid for the 
free electron case, but it is in error for the non-spherical Fermi surface 
due to the contribution of some terms which the adiabatic approximation 
necessarily neglects. However, as we shall restrict ourselves mostly to 
the free electron case, the adiabatic approach is still a sound method of 
calculation in what follows. 
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We assume that the electron-electron interaction, expressed in the 


interaction representation, is exponentially switched on between t= — © 
and t=0, i.e. is represented by 
H(t) =exp(it t)H, exp(—tH ot) exp(t) . = (2.14) 


where the parameter « governs the rate of switching on the perturbation. 
The adiabatic theorem (Gell-Mann and Low 1951) establishes that if 
ly.,(¢)) satisfies the Schroedinger equation 


iM IQ) = 20D) + + RNB) 


and initially (t= — 00) is an eigenstate of Wp, Le. 


lx..(— 20)) =|), 
and if in addition the limit 
J (ELA a |x.(9)) é 
AF= lin —_—“— Pe as See 
a0 (‘E 1x2(0)) 
exists and is real, then the state-vector 


WY — fim X29) 

Ley lim HT y.(0)) ee ee OES ye 
is an eigenstate of the total hamiltonian # with an energy eigenvalue 

E=W-+AE, 

where W is the unperturbed energy given by (2.13). This means that 
when the interaction is switched on infinitely slowly, an eigenvector of the 
unperturbed system is transformed into an eigenvector of the total 
hamiltonian. 

Assuming the existence of the limit (2.16), if the unperturbed state 
|W’) is non-degenerate, AH is always real and the hypotheses of the 
theorem are satisfied; this is the case of the ground state |@). Lf, on the 
other hand we apply (2.16), (2.17) to a state |’) which is degenerate, then 
AE turns out in general to have an imaginary part, its physical meaning 
being that the perturbation transforms the original eigenstate |’) into a 
decaying state |’) of the total system with a mean life given by 

T={lonART: ogee! se eet) 
After a time of the order 7', the system decays into other states due to 
electron—electron collisions. It is theoretically possible to find linear 
combinations of degenerate states |’) such that the imaginary part of 
AE vanishes for each of them; these are the true excited eigenstates of 
the total hamiltonian. However, due to the formidable degeneracy of 
the electron gas and the enormous difficulties in exactly evaluating the 
matrix elements (2.16), the problem of finding the exact excited eigen - 
states has never been solved. Even if we could calculate them explicitly, 
these would probably be of very little use for studying the transport 
properties of the system ; the complicated linear combinations of different 
|") would mix the states in such a disordered way that no one-to-one 
correspondence between the eigenstates of the independent particle 
model and the true hamiltonian would be possible. 
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However, as discussed in § 2.3, it is possible to find some unperturbed 
excited states which give rise to decaying functions |’) with a very long 
mean life, i.e. with a very small imaginary part of AH. If one such state 
is excited by an external agent and its mean life is longer than the relaxa- 
tion time of the whole process (due to impurities or the motion of the 
lattice), it can be considered for all purposes of transport theory an eigen - 
state with energy 

H=W=ReAE. 


The imaginary part of # can thus be taken into account as a small correction 
to the relaxation time. 

To compute the energy shift (2.16), the operator S,(t) (similar to the 
S-matrix of quantum field theory) is introduced in such a way that 


(EVI) =) eat eee ee, eo (219) 
It satisfies the Schroedinger equation 
iS MLE oe ALS AG oem ot Pe) 8s fin 8s = ye (2020) 
with the initial condition 
S,(— c)=1. 


Solving (2.20) by iteration under the assumption that the perturbation 
series converges, the final formula for S,(0)=S, is obtained : 


Seige Ss (aye ml wide! ee (ae it, 
n=1 
x PL H(t) HL, (ty)... #,( erro t (2/91) 


where the operator P places the eee ¢ in Pi enclsaien order such 
that the times never decrease going from left to right. 

Inserting (2.19) in (2.16) and (2.17), we write the perturbed energy and 
state vector as 
(FLA4S.1P) 


(EIS, |) 


a0 


and | 
; Si|'¥) 
aie CaS, Wg 

Finally substituting (2.21) for S, gives perturbation series for AH and 
|’). Each term of these series can be represented by a diagram similar 
to those introduced by Feynman (1949) in quantum field theory. For 
the sake of definiteness, we give explicitly the prescription for building 
up one of the graphs contributing to the numerator of (2.22), and the rules 
for calculating from it the corresponding contribution to the energy: 

(i) Mark n+1 pairs of points on the diagram ; both points of each pair 
must be at the same time (considered to increase from bottom to top of 
the diagram), but different pairs are at different times ; 

(ii) Label each pair with times fo, t,,...t,, beginning with the uppermost 
pair and proceeding downwards ; 

(iii) Join both points of each pair with an interaction (dashed) line ; 
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(iv) Draw particle (full) lines joining different points, or going from one 
point to the upper or lower end of the diagram, or going vertically across 
the graph without touching any point; put arrows on them such that at 
each point there are two and only two particle lines, one with incoming 
and the other with outgoing arrows, and at both ends of the diagram 
there are as many lines pointing upwards as electron-like quasi-particles. 
exists in |’) and as many lines with arrows pointing downwards as hole- 
like quasi-particles exist in |’). 

(v) Label each external line (line going to one or both ends of the 
diagram) with indices k, or k, corresponding respectively to the electron- 
like (upwards arrows) or hole-like (downward arrows) quasi-particles of 

pele 
| (vi) Label each internal line (line joining two points of the diagram ) 
with indices k, k’,k”, etc. 

To evaluate the contribution of each diagram to the series for the 
numerator of (2.22): 

(A) Write for each interaction line the factor 

20234 EXP[ZU(ky? + ky? — ky” — ky”) jt, exp(at,,) 
where k, and k, are respectively the labels of the outgoing and incoming 
particles at one point of the pair ¢,,, and k, and k, the labels of the out- 
going and incoming particles at the second point of the pair; 

(B) Write a factor 

Dt, ar tm) = Y(t, - tm) =F Y(k ay ky) 
for each internal line k going from a point ¢,, to a point ¢,; the function 
Y(x) is defined by 
Yit\jent axe G, 
=I) Facets: 

(C) Put t)=0 and integrate over all other time variables between 
— o and 0; 

(D) Multiply the integral by a factor 

(=i"(n (= 1 
where n is the number of interaction lines minus one and y the number of 
closed particle loops of the diagram ; 

(EK) Sum over all k, k’, k”..., the labels of the internal particle lines ; 

(Ff) Evaluate the limit «+0. 

Similar prescriptions can be given for all the other relevant matrix 
elements. The results which follow in the next subsections have been 
obtained by systematic use of this or other very similar techniques. 


2.2. The Ground State 
In evaluating 
, 4 GIFS 1G 
AE, = lim (Gl eeeetG) of ean te Yee A ee Ae 
x->() (GS,|G) 


Goldstone (1957) has proved that disconnected diagrams such as that. 
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shown in fig. 3, i.e. diagrams with two or more different parts not connected 
to each other by lines, give contributions only to the normalization of 
the wave function. Although the integrals arising from these disconnected 
diagrams diverge, their contribution is exactly cancelled by the denomina- 
tor of (2.24); they can be completely disregarded and the series is then 
reduced to only connected diagrams. This is called the linked-cluster 
expansion. The linked-cluster perturbation formula was also proved by 
Goldstone to be a continuous function of « and therefore the limit can 
be taken by putting w=0. The final formula for the energy is expressed by 
1 

Alem BC yy, 
where >, means summation over connected graphs only. The same result 
was obtained by Hugenholtz (1957, 1958) by a somewhat different method. 
To evaluate (2.25) explicitly it was necessary to overcome a crucial difficulty: 
the divergence of the summation over the intermediate states k, k’,... 
due to the piling up of factors 1/k?. The first-order term of the series is 


#;)"|6) Poh ees25) 


Fig. 3 Fig. 4 


2 


shown in fig. 4 and its calculation is straightforward ; it gives the exchange 
energy in the Hartree-Fock approximation. The second-order diagrams 
are shown in fig. 5; the first one is the direct term and diverges logarith- 
mically, while the second represents the exchange term and is finite. 
Gell-Mann and Brueckner (1957) have proved that the divergence of 
the direct term is cancelled by similar divergences in high-order terms 
such as those shown in fig. 6. By summing the contribution of these 
diagrams up to infinite order before performing the summation over the 
divergent internal line they were able to get an expression for the energy 
of the ground state in terms of the interelectronic spacing : 
EL, = N[2:21rs — 0-916 rs + 0-0622 In rs — 0-096 + O(rs)Iry. . (2.26) 
This expression is only valid in the high density limit, rs <1 as has been 
shown by Noziéres and Pines (1958), and is therefore not useful for metallic 
densities 2<7rg < 5:5. 
Hubbard (1958a) using the effective interaction method developed 
by him (Hubbard 1957) was able to repeat Gell-Mann and Brueckner’s 
calculation (to second order in the perturbation expansion). He also made 
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an approximate evaluation of the third-order diagrams, obtaining results 
which may be valid even in the region of metallic densities. These are 
shown in table 1, where values of the formula derived by Bohm and 
Pines (Pines 1955) using the collective description are also displaced. 
The latter is given by 

By = N[2-21rs2—0-916rs-1 + 0-031 Inrs—0-115]ry . . (2.27) 
and can be considered a suitable interpolation formula for metallic 
densities (Noziéres and Pines 1958), where the last two terms give the 
correlation energy within an accuracy of about 15%. 


Fig. 5 


Table 1. The energy of the ground state 


Ke 


Hubbard’s value Bohm and Pines’s value 


— 0-099 — 0-093 
— 0.086 —0-081 
— 0-074 — 0-072 
— 0-067 — 0-066 


Note : Eg is usually called the correlation energy and is defined by 
Hg=N[2-21 rs —0-916 rg? + Eo]ry. 
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Concerning the structure of the ground state, i.e. the explicit form of 
the eigenvector |G), the results obtained are very few and only related to 
some particular details. The reason for this partial failure of the method 
lies in the fact that the series (2.23) for the state vector, due to the divergence 
of the normalization terms, is essentially non-convergent ; the cancellation 
of terms that give rise to the linked-cluster expansion for the energy is 
not valid here and up till now every perturbation method has failed in 
giving even an approximate solution for |G) (see, for instance, Brueckner 
1958). However, Migdal (1957) and Luttinger (1960) have proved that 
the average occupation probability f(x) in the ground state |G), although 
different from the simple step function characteristic of the unperturbed 
|G), still exhibits a discontinuity Af<1 at some values of k, in terms of 
which the Fermi surface can be clearly defined (fig. 7). Luttinger (1960) 


Fig. 7 
(hk) f (R) 


M4 R 
~k, 0) Re F 


(a) (b) 
The occupation probability f() for the (a) unperturbed and (6) perturbed 
ground state, showing discontinuity at the Fermi level. 


has also proved that the volume in k-space enclosed by the Fermi surface 
defined in this way is equal to the volume of the Fermi distribution of the 
corresponding non-interacting system, i.e. contains exactly 1 or 2 or... 
electrons per atom. 


2.3. The Single Quasi-particle States 
Suppose that in (2.22), (2.23) we put |’) =|k,), i.e. a single unperturbed 

electron-like quasi-particle in an N-+1 particle system (or a hole-like 
quasi-particle in an N—1 particle system). Due to the degeneracy of 
unperturbed states having this total k and total energy, the energy (2.16) 
turns out in general to have an imaginary part as already mentioned and 
the adiabatic theorem that (2.23) is an eigenstate collapses. Nevertheless 
we can still use (2.23) to define the perturbed state |k;) corresponding to 
|k;) cats 

k= lim 7% 6 www we (2.28) 

| ? a>0 (k,|S,,|k;) 
and investigate its properties, 


78 L. M. Falicov and V. Heine on the 


Consider first the case k,=k, with the quasi-particle just at the Fermi 
level. The degeneracy of the unperturbed state tends to zero, and Hugen- 
holtz (1957, 1958) has proved that the adiabatic theorem holds and we 
obtain an eigenstate with a real energy 

E(k) = Egy + (kK) Gos eee 
where the plus or minus refers to electron-like or hole-like quasi-particles. 
He has also shown that the state |k;) defined by (2.28) can be written in 
the form of an excitation of the ground state, 


|k,, V-+1)=OMIG, NV). 8a eel 2-e0 ay 
where 
=2 n 
oh fq, \* a lee ae Mek 
O, Ny 2,2, Ut, s boa -) | i ( ) 


Here {a,'* is the product of creation and annihilation operators of the 
‘bare’ electrons that transforms |G,.V) into an arbitrary unperturbed 
state |v), and >, means that the sum over n must be performed taking 
account only of connected diagrams whose intermediate states are all 
different from |k,). The energies in the denominator are the perturbed 
energies. Similarly for holes 

Ik. V1 SOIGSN). 4). oer ae 


where O,, is the hermitian conjugate of O,*. 

For k,>k, (or k,<k,), we can use either (2.28) or (2.30) formally to 
define the perturbed states. Incidentally it has only been proved that 
the two states so defined are identical at k,: away from k, there may be 
slight differences in the imaginary parts increasing with k,—k,, but for 
practical purposes near the Fermi level we shall treat them as identical. 
We shall use (2.30) formally to define our states as in § 1.2. but revert to 
the series for (2.22), (2.28) for quantitative calculations. The energy «(k), 
still defined by (2.29), varies as follows (Hugenholtz and van Hove 1958): 

Re e(k) — €, ck—ky, \ ee ae 
Ime(k) c(h —ky)? 
where we define ¢, as the Fermi energy. Hugenholtz (1957, 1958) showed 
that in evaluating e(k) only connected diagrams need be taken into account : 
the disconnected diagrams contribute only to E,, y and to the normalization 
of the state vector. 
The Fermi energy ¢, has the property 
Bes vi &y= Hg, x =Heg vite 
showing that ¢, can be defined as the separation energy of one electron. 


In fact this equation as it stands is not very meaningful because Eg y is 
of order N and N->+oo. It can be rewritten in the form 


e 2 [Xe 
Pwdp\ Q 


which has a well-defined meaning for an infinitely large system. As usual 
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the Fermi surface is defined by ¢(k)= e, and the whole calculation, 
particularly (2.32), shows that its radius &,, is still given by (2.10) the same 
as in the Sommerfeld model. This shows it is also the same as the Fermi 
surface defined in §2.2 by the discontinuity in the occupation function 
so that no ambiguity can arise. 

The operators O,*, O, are creation and annihilation operators for 
‘dressed’ quasi-particles. Thus as described in §1.2, we can set up a 
whole system of excited states of the electron gas in terms of multiple 
excitations of the ground state. From their mode of derivation, they 
must at the Fermi level satisfy the usual anticommutation relations for 
fermion operators, although their structure is too complicated to verify 
this explicitly. It is not clear whether the same property holds when the 
definition (2.31) is extended away from the Fermi surface, although it is 
probable that this is not so, with the consequences discussed in § 1.3.3. 

Quinn and Ferrell (1958) have computed ¢(k) to second order of pertur- 
bation theory using the effective electron—electron interaction given by 
Lindhard (1954) and eliminating the divergences by subtracting the 
contribution of the vacuum. Although the calculation is only valid for 
very high densities, it shows some important features of «(k) and the 
lifetime of the quasi-particles. As far as we know, no good calculation of 
e(k) for metallic densities has been completed up till now. In fact all the 
available formulae are only rather crude approximations. In this paper 
we shall take for performing explicit calculations, the Bohm and Pines 
formula (Pines 1955) obtained through a quite different approach, i.e. the 
collective coordinates model. For quasi-particles near the Fermi surface 


0-611 3k? + (B2—1)kep? 
ibs a es es ee ees 2 eee 
c(k)=| ar {1 28 + a 
ky?—h? kp +h 
ewes ree oes ee 33 


where f is a parameter that governs the cut-off of the effective interaction 
(see §2.4). According to Pines (1955) 
Bassas eee) Dako 0 L984) 
and Fletcher and Larson (1958) give 
B=0-4r,1? 
as a better value that brings the theoretical values to a better agreement 
with experimental results. 

A more accurate calculation of ¢(k) would be desirable and even necessary 
to apply theoretical results to the real metals ; but the discouraging mathe- 
matical difficulties of the problem have restricted the application of well 
elaborated methods, such as Hubbard’s perturbation techniques, to 
numerical evaluation of the ground state energy only. 


2.4. The Effective Quasi-particle Interaction 


From the work of Bohm and Pines (Pines 1955) it is well known that 
the collective properties of the electron gas result in an effective screening 


P.M.S, G 
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of the coulomb interaction. It has been shown that the long-range part 
of the interaction results in some collective oscillations of the system as a 
whole; these collective oscillations have a very high frequency and are 
excited only in very energetic processes. The remaining short-range part 
constitutes the effective quasi-particle interaction. Bohm and Pines, in 
a somewhat arbitrary way, eliminate the long-range effects by changing the 
Fourier transform of the coulomb potential 


v(k) = 42/k? 
into an interaction 
=0, k< Bk, 
Vpp(k) 4 Ply ie ae Pe Dean 
| =4n/k?, k>Bky 


that shows a cut-off on the small k-vector side. The cut-off parameter fh 
was chosen by a variational argument, but many doubts have arisen 
concerning what is the most suitable magnitude to minimize ; the selection 
of the best 6 is not yet quite clear (see for instance Fletcher and Larson 
1958). In a way it is a rather academic question anyway since the form 
of the interaction is too unrealistic for most quantitative calculations (see 
below and §5). 
Fig. 8 
¥--- ¥---- 


---X ----x 

The screened character for the effective interaction and the existence of 
collective motions of the system have been also obtained by means of the 
many-body approach. Hubbard (1957) has shown that some parts of the 
diagrams in the perturbation series can be summed up to give a modified 
interaction, which is time-dependent and short ranged. 

The portions of the diagrams to be summed are the so-called polarization 
parts and consist of those sections of any diagram connected to the rest 
only through two interaction lines. The summation, that is schematically 
shown in fig. 8, means that all polarization parts of the diagrams can be 
ere if the interaction v(k) is replaced by an effective potential V. 

ut now, since the interaction lines join points at different times. V S 
necessarily be time dependent, yeehanta® a retarded meine "The 
Fourier transform (f to w) of the effective potential can be expressed ast 

v(k) 
e(k, w) 
where e«(k,w) can be interpreted as a dielectric constant. energy and 


V(k, w)= (2.36) 


| The double use of ¢ for one-particle energy ¢(k) and for dielectric constant 
e(kw) in this section only, should lead to no confusion. 
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momentum dependent, due to the electron gas. Hubbard (1958 a) has 
computed <«(k, w) up to second order (in first order «=1) and his result 
agrees with the formula given by Lindhard (1954) and used by Quinn 
and Ferrell (1958) to compute single-quasi-particle energies. ‘This second- 
order dielectric constant is a good approximation for the high-density 
limit, but not suitable for real electron gases in metals. None the less it 
clearly shows the screened character of the interaction. Putting 


eK ch) mel I KS a)! aaa) ot iat > 2 (2:37) 


the second-order computation gives 


k ( w nite ) i 
F OTe) yee > co cap Ty ED 
F i(k, w) = un A ee) 2hy ] w — kky — 3h? 


mk? k w+kk, — tk? 
w k \2 
be (Ge + 3p) 6 is wtkkpthk2 
cs i Se 


This, for w=0 and k very small, reduces approximately to 
Fi (k, 0) — 4hy/ak® i ci at ee (2 80c) 


which is equivalent to an interaction V(k) equal to a Yukawa type of 
potential (1.8) with effective screening length 


1/Q=(a/4ky)¥2. (2.39 6) 


The approximate expression (2.39) has been used by Quinn and Ferrell 
(1959) to study some properties of an idealized metal at high density and a 
similar expression has been used in § 1.3 as prototype of an effective screened 
interaction. 

We must also note that for small & and large w, (2.38) reduces to 


F,(k, w) » 4k,3/37w? 


showing the possibility of collective oscillations of the system with the 
classical plasma frequency 
wo, = (4h,,2/3)"? 


for which the dielectric constant vanishes. 
For metallic densities, it is necessary to evaluate ¢(k, w) at least up 
to third order of perturbation theory. Hubbard (1958 a) has evaluated 
- approximately the contribution of the third polarization sub-diagram 
shown in fig. 8 and a series of related diagrams, with the result 


F4(k, w) 


[Wee 
ee AK 
45 e+ key?) al w) 


F(k, w) = 
1 


(2.40 a) 


TG 2 
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A better formula is probably (Hubbard, private communication) 
FA(k, w) 

L+ 3k? F'o(k, w)/(k? + kp? + Q") 

where Q is given by (2.396), though neither of these has quite the correct 
form for small k& (see fig. 9) where (2.38) is more correct. The different 
results are plotted in fig. 9, where we also include our fudged best 
estimate based on (2.39) and (2.406). All of these calculations neglect 
certain effects which must be of some importance at actual metallic 
densities, so that one cannot say how reliable the results are. 


F(k, w)= (2.406) 


Fig. 9 
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Various approximations for V(k,0) for r,=3 using: (a) Bohm and Pines’ 
cut-off function (2.35), (2.34) ; (b) Lindhard’s second-order calculation 
(2.38) ; (c) Yukawa interaction (1.8), (2.39); (d) Hubbard’s original 
calculation (2.40 a). The fact that the curve tends very nearly to zero 
at k=0 is not a general feature but is associated with the particular 
value r,=3; (e) our best estimate (after Hubbard, private communica- 
tion). Within the limit of the effects being considered here, this should 
tend to Lindhard’s value at k=0. For k21 we have used (2.40 b) 
(shown dotted for small &). The ordinates are in units of Ar |Iop?, 
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The introduction of the quasi-particle formalism and the effective 
interaction, allows us to make a change of representation of our hamiltonian, 
which now becomes 

WO ee en ent) ar abides. oo! + (2.41) 
where 
a = Ky aii 2, k;)O,;* ‘Og + > ¢( On,O, z 


oh 
gives the single quasi-particle energies, and 
eas 2 Faaat Ora “O23 “Orca} ev eto) 
gives the quasi-particle interaction term. Here 


V iosa oT (1/Q)V(k, = k,, Seon €4)0(4, 5 4)0(Oy, o3)0(ky +.k, ks aD k,) (2.44) 


Fig. 10 


and >* means that in calculating the evolution of the system, diagrams 
with interactions not connecting different quasi-particles as well as those 
containing polarization parts, must be disregarded. The bracket means 
that in ordering the O operators, the pair O,,*O,,, for holes must be replaced 
by —O,,0,.*, because the contribution coming from the term | in the 
commutation relation is already taken into account in the ground state 
energy Hy. 

Since we are only dealing with quasi-particles near the Fermi level, 
we can put e«,—«,=0 and consider only the so-called static effective 
interaction V(k,0). Figure 9 shows various approximations for V(k, 0) 
for a typical metallic density of r,=3. We note that our best estimate 
(curve e) is much smoother and considerably smaller in average magnitude 
than the Bohm and Pines interaction. We shall return to this in $5 
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when discussing numerical estimates of the correction terms due to the 
interaction. 

Consider now a state containing two quasi-particles, The interaction 
H ,' (2.43) contributes in first order of approximation an interaction 
energy 

71 (Ky, Ky) = (Vi201 + Viera) 

which amounts to + (1/Q)V(k,—ky,, 0) for particles of parallel spin and 
zero for antiparallel spin. The minus refers to two electron-like quasi- 
particles or two holes, the plus to an electron interacting with a hole. 
We can now apply %,,1 in higher order of perturbation, represented by 
diagrams of the sort shown in fig. 10, and the second order ones give small 
contributions 7, both for parallel and antiparallel spin. The total inter- 
action energy 7(k,, ky) of (1.13) is taken to include the effect of all orders 
of perturbation in this way. 


2.5. The Effect of the Periodic Potential 
To find the eigenfunctions and eigenvalues of (2.1) when 
U(G) 40 


implies the solution of an extremely cumbersome problem that has never 
been solved up to the present, even for simplified or trivial models of 
lattice potentials. However, the existence of the Fermi surface, defined 
as the points in k-space where occupation density f(k) is discontinuous, 
has been proved very generally by Luttinger (1960) and this, as well as 
the generality of many conclusions drawn from the free electron case, 
seem to indicate the existence of quasi-particles of long mean life near the 
Fermi surface. 

Of course, the existence of interband transitions in the metallic band 
case can bring important quantitative correction to the properties of the 
quasi-particles and their interactions. However it is reasonable to assume 
that their mean life decreases quadratically as they approach the Fermi 
level while their energy varies linearly, which is all that is necessary for 
the quasi-particle picture to hold. It is also very probable that the effective 
interaction, although it must exhibit some angular dependence, is not 
too different from the free electron value in the sense that it is very unlikely 
to exhibit strong oscillations or sharp differences, because it really is an 
integrated effect involving all the electrons. Of course in metals with 
incomplete d shells the effects of correlation can become very large and 
even lead to a Heitler—-London localized type of wave function. 

All these assumptions must be confirmed or corrected when definite and 
complete calculations are available. Up till now only methods of calcula- 
tion have been proposed, but as far as we know none of these has yet been 
carried through explicitly. Hubbard (1958 b) has developed a modification 
of his free electron treatment to study the ground state. The method 
involves a formidable amount of numerical computation to solve, firstly, 
the Hartree equations for the electrons in a lattice, and secondly the 
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effective interaction, expressed as an integral equation, when the screening 
is modified by the ionic field. The method has been extended by Pratt 
(1960) to obtain a Schroedinger-like equation for the quasi-particle 
energies ¢(k). 

The formal difficulties noted by Kohn and Luttinger (1960) and Luttinger 
and Ward (1960) about the applicability of the Goldstone formula have 
already been mentioned in § 2.1. 


2.6. Derived Properties 


Quinn and Ferrell (1958, 1959) have studied various physical properties 
of an idealized metal, using a quasi-particle approach similar to that used 
in the present work. Their results are qualitative and applicable to the 
high density limit only. 

Wolff (1959) has studied the effect of the electron correlation on some 
optical properties of metals. He has used the diagrammatic analysis of 
Hubbard and has obtained a formal expression for the current when the 
electron gas is excited by an external electromagnetic field. In the free 
electron case, no contribution from the electron—electron interaction 
appears in the current as expected. However, for non-spherical Fermi 
surfaces the contribution can be important and Wolff discusses the metallic 
gas in the extreme and idealized case of two isotropic bands of electrons 
and holes respectively. 

The thermodynamical properties of the electron gas have been the 
subject of several contributions. In particular Gell-Mann (1957) has 
calculated the specific heat of the free electron gas in the high-density 
limit by a method similar to the one which led to the elimination of the 
divergences in the ground state energy due to the coulomb interaction. The 
determination of the thermodynamical functions of the interacting electron 
gas by means of diagrammatic analysis has been published by Lee and 
Yang (1960), Montroll and Ward (1958) and Bloch and de Dominicis 
(1958). The technique used by the latter has been employed by Luttinger 
(1960), Kohn and Luttinger (1960), and Luttinger and Ward (1960) to 
test the validity of Goldstone expansion in an arbitrary non-free electron 
case, as well as to calculate the specific heat and other properties. 
Brueckner and Sawada (1958) have calculated the paramagnetic spin 
susceptibility in the high-density limit. The specific heat has also been 
calculated by Chen and Chow (1958). 


§ 3. VELOCITY, CHARGE AND ACCELERATION 


When we ask what is the current carried by a quasi-particle an apparent 
paradox immediately appears. One answer would seem to be 


ik Steve eee pot eel eS Ac) 
e(iK ‘ 
where v(k) = ; ee ) (3.2) 


is the usual group velocity of the quasi-particle and ¢ the electronic 
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charge. On the other hand, from fundamental quantum mechanical 
arguments 


jm (ey, eae eee ec 


m 
which relates the current to the expectation value of the total momentum 
P. In the case of a single electron, free or in a periodic potential, it is 
easy to prove that both answers are the same (Wilson 1953, eqn. (2.8.3.)). 
However, we can see that (3.1) and (3.3) are not equivalent in a many- 
particle system when the mutual interaction is taken into account. The 
discrepancy is evident even in the simple Hartree-Fock aproximation 
for the free electron gas. In this case the momentum of a single particle 
is Ak, but 
e(k) A h7k?/2m 
due to the contribution of the exchange term, and therefore 
ev Zenhk/m. 
We shall resolve this apparent paradox in §3.2 through a calculation of 
the effective charge e* carried by a quasi-particle. 


3.1. The Velocity 


The usual argument of the group velocity of a wave leads to the well- 
known formula (3.2) for the velocity of a single quasi-particle. However, 
since the definition of effective charge (§ 3.2) rests mainly on this concept, 
it is worth giving a detailed argument. 

We shall need the single quasi-particle states |k;) of the V + 1 electron 
system defined in § 2.3, and their development in time 

|k;, ¢)=|k;)exp[—w(Hyte,)t/h] . . . . . (8.4) 
where we have assumed that ¢ is sufficiently long to be of physical interest 
but short compared with the natural lifetime 7’ of the state so that we 
can treat |k;) effectively as an eigenstate of the hamiltonian. 


To study the evolution of the system in real space, we introduce the 
density -of-particles operator 


p?(r) =(1/Q) ¥ a,*a,.4expiq:r. «ows Ae a ae 
k,q 
The mean value of this operator for a state |Y’(t)) 
p(r, t) = CF (t)|pP(r)[F(t) > 
gives the density of particles at r at time ¢t. We observe that in (3.5) each 
term only connects states with total momenta differing by an amount q. 


Therefore, the terms contributing to ¢k;|p°?(r)|k;) are only those with 
q=0, and the particle density for the state |k;) 


p(r, t)=(1/Q)¢k;| Sa,*a,]k;) = (NV + 1)/Q 
is uniform. This is, of course, a trivial result since the electrons must 


necessarily be uniformly distributed when the total momentum is a good 
quantum number. 
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In order to get a localization of particles in real space, we form a wave 
packet of quasi-particles 


|ko, a) = >,a(k;)|k; > . . . . . . . (3.6) 
where a(k;) #0 only in a small region of k-space around k, and 
Diep) [eae Lowen el -arey 36 deere Ait is Soo (3875 


so that the packet is normalized. 

We can always choose the coefficients o(k;) in order to obtain a high 
density of electrons concentrated in a small region about ry at a given time 
ty. The density of electrons for this wave seoihe is 


p(r, t)= (Ko, a, tp? (r)|ko, «, £) = 2oM(k 5) Pri, jist) ©. (3.8) 
Pres, wiih #) = ¢k;, t[p°?(r) i, b> Riley, ona. tot(3.9) 
These matrix elements can be written in the form, 
PARE) Otel La pits. ne Sa. (S10) 
Pri, ki-+-q = B(k,, Renee es on, +q)—e(k;)}é/h, q#0, (3.11) 
Bik; k,+q)=(1/2)>,<kila,*a,, glk; +4)... . (3.12) 
Inserting (3.10) and (3.11) ah ae we oe 
N+1 
pr, )=—2— + Ya (k,)a(k; +4) Blk; ki +4) 
Or ki, q 40 


: l/de aay 
xexp| iq.{r—7(55), Seen (cbalcn) 


where, using the fact that a(k,) is different from zero only around k;= kg, 
we have made the approximation 


cl) = (ko) + (hie) (5) 


(3.13) shows clearly that the high concentration of electrons, which at 
t=t, was centred at r=ro, at a later time ¢=t,+ At is centred at 


r=r,t+Ar=ry+ + (5) At 


and consequently is travelling with a velocity 


Ar l/de 
V(Ko) = ee (i), 


as we wanted to prove. 

When we take account of the lattice potential we label our state |k,, s) 
by the reduced wave vector k; and band index s. Each term of the particle 
density operator (3.5) now connects states whose k-vectors differ by q, 
where 

k’-k=q+G, 
and G, is any reciprocal vector of the lattice. Forming as before the 


wave packet 
ko, 8, a) = Da(k,)[k;, 5) 
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the mean value of the density of electrons for this packet is still defined 
by (3.8) and (3.9), but the matrix elements are now 
Pug ee eG Kor): . 2 eee gh ocea) 
Pri. ki = Uk; ki +4; r)expliq. r—if{e(k;+q)—e(k,)}t/h], (3.15) 
where the functions M are periodic functions of r with s periodicity of 
the lattice. For phenomena involving many cells of the crystal (and this 
is the case we consider throughout this paper), the spatial modulation 
can be ignored by replacing M(r) by its mean value M. It follows, then, 
that 
M(k;, k,) =(N + 1)/Q, APE ar A GRY 
M(k;, k;+4q)= ote >, (ky, 214,"4,.., Kendo), eee 
and the argument given above for the free electron case can be applied 
directly. 


3.2. Current and Effective Charge 
From fundamental quantum mechanical theory, we know that the 
current carried by a system of many particles is given by the expectation 
value of the operator 
Se 
7 ™,; 


z 
which in the case of the electron gas is 


where P is the total momentum operator. Thus in the case of a single 
quasi-particle the current is defined by 


j(k) =(k | 


m 


k). foi & ee ees) 


As already mentioned, we cannot expect j(k) and ev(k) to be equal. 
It is therefore convenient to introduce an effective charge e*(k) for the 
quasi-particle defined by 


ark) = MW CGA er ee ERE 


so that the charge e* travelling with velocity v gives the correct current j. 
We shall prove in the appendix, inter alia, that v and j are always parallel 
so that e* is a scalar. The concept of an effective charge has also recently 
been developed by Stern (1960). 

We shall now give a physical interpretation of e* and show that if we 
form a wave packet, the actual number of electrons carried along by it is 
in fact e*/e. Suppose we have a packet of SN electrons localized 3 in gue 
at a particular instant of time, say (=0. The particle density p(r) is 
then as shown in fig. 11 (a). If we express the particle density as a Fourier 
series 


p(t) = >,4(q) exp(iq:r), 
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then A(q) is a continuous function of q at q=0 and the number of electrons 
5N is given by 
ONS QA(O)Sa lime OA(q)e suit 5. (3.20) 
q>0) 
Suppose now we add a constant uniform electron density (V + 1—8N)/Q. 
Then p(r) becomes (fig. 115) 


p(r) = ages eq) expiiqur)— . =. (3.21) 
N+] 
=- 5 tay eA (Q)OxT eG 0 ).c in ine ce) 03.22) 
Fig. 11 
e 
5N (Qa) 
0 ete 
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Particle density in a wave packet. 


where >’ excludes the term q=0. It is now clear what we must do to 
calculate the number of electrons 6N transported by a quasi-particle, 
We already have in (3.13) the charge density for a one-quasi-particle 
wave packet, expressed (at any particular instant) in the form (3.22). 
In order to calculate the number of electrons 5N localized in the bump of 
fig. 11(b), we have to work back to the form (3.21) for p(r), and thus we 
obtain i 
SN = lim QY,,0°*(k,)o(k;+q)B(k;, kj+q). . . . (3-23) 
q>0 
In the Apprendix this expression is evaluated and it is shown that it gives 


SN =e*(k,)/e 
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where e* is defined by (3.19). Now as the wave packet travels with time 
(fig. 11 ¢), the charge transported is e 5.N = e* and the current is) =¢dNv= e*v 
which is our interpretation of (3.19). In order to complete this interpreta - 
tion of (3.19), we have to show in addition that the uniform background 
p(r)=(N +1—5N)/Q is in fact stationary and does not contribute to the 
current; this is also done in the Appendix. 

In the free electron case, the perturbed state |k) is still an eigenstate 
of the total momentum P since P commutes with the coulomb interaction, 
and the current from (3.18) is 


jk) = chk fics Si eee 
The effective charge is therefore given by 
Rk FR | Cea 
Me oe AT. . . . . = = ee 
cull mv mM (5) re 


Table 2 gives some values of e*(k,)/e in the Bohm and Pines approxima- 
tion, and we note they are considerably less than one. We expect a more 
realistic interaction to give values much nearer unity as discussed in § 5. 


Table 2. The effective charge of quasi-particles at the Fermi level 


2 e*(kp)/e 
0-2 . | 0-98 (a) 
0-5 0-98 (a) 
2-0 0-87 (b) 0-90 (ce) 
3-0 0-88 (b) 0-92 (ce) 
4-0 0-90 (b) 0-95 (ce) 
5-0 0-93 (b) 0-96 (ce) 


(a) High-density limit. 
(b) Bohm and Pines’s approximation with 8 =0-353r 12. 
(c) Bohm and Pines’s approximation with 8 =0-47 1. 


3.3. Acceleration in Applied Fields 
A single quasi-particle, which we may picture in terms of a wave 
packet, moving in an electric field & picks up an amount of energy 
de=e*&-dr=6: jdt 
where j is the current carried by it. Thus the acceleration in k-space is 
given by 
ree de/de e*&-v eG 
=—/~ =———_ =~ .... . (8,26) 
dt/ ok hv h 
In practice of course one always has multiple excitations, and it is much 
more difficult to calculate the acceleration then because some of the energy 
picked up by one quasi-particle goes into the interaction energy with the 
others. It is for this reason that we avoid any such attempt to calculate 
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the total acceleration in our derivation of transport coefficients, and use 
the method of § 1.3.4. We may use (3.26) as an aid to the imagination in 
guessing the type of multiple excitation relevant to a conduction process 
(not its magnitude) and this may be somewhat in error, but our method 
(1.20) of evaluating the resistivity from it is very closely related to the 
variational principle of transport theory as discussed in §5, so that the 
errors in the answers are expected to be quite small. 

By analogy with (3.26), we assume that the acceleration in a magnetic 
field is given by 


ik= (iclevaHe o 6 3 ok 2 38n 


Since the derivation of (3.26) was based on the conservation of energy, the 
same technique cannot be applied to the magnetic case. However, the 
e* seems to be demanded in (3.27) by the macroscopic law of force ja H/c. 
Also by the relativistic connection between electric and magnetic fields, 
we expect the same factor in (3.26) and (3.27) ; at least it would be difficult 
to justify the use of e instead of e* in (3.27). We have set down these 
rather obvious remarks, because there is considerable uncertainty about 
the magnetic effects (§§ 4.6 and 5). 


§ 4. ‘ TRANSPORT ’ PHENOMENA 


We can now discuss those properties of a metal which depend only on 
excitations at the Fermi level. These include the electronic specific heat, 
paramagnetic spin susceptibility and de Haas—van Alphen effect, as well 
as most genuine transport processes. As already pointed out in § 1.3.4 the 
latter must be treated by the method of energy dissipation. Now most 
transport phenomena involve in an essential way the collision processes 
suffered by the electrons, and the general uncertainty about these is so 
large that it completely overshadows the kind of few per cent corrections 
which we are now discussing. Thus we shall restrict ourselves mostly to 
those effects which are directly related to the shape of the Fermi surface 
and the Fermi velocity, independent of the realxation time rv. 


4.1. Electronic Specific Heat 


Calculating the electronic specific heat is fundamentally a problem in 
statistical mechanics, for which we have to enumerate all the low-lying 
states of the system with their total energy, so that we have in principle 
to consider the energy corrections of §1.3.4. Now at a temperature 7’ 
we have of the order of N&T7'/e, excitations, i.e. a number of order J. 
In an overwhelming porportion of states, these excitations are more or 
less uniformly distributed round the Fermi surface with fluctuations of 
order \/N. The total interaction energy (see eqn. (1.15)) 


H,= >> (ky, kg) 


is then zero because for every 7(k,, kz) between two electrons or two holes 
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there will be an 7 of opposite sign from an electron and a hole with practically 
the same k,, k,. Only the individual particle energy 

Ey = Ye(k;) — De(k,) 


contributes, and the electronic specific heat is the same as that given by 
the independent electron model, namely 


O,syL= tla <a): eee ee 
1 dS 1 dS 
‘ N (ep) = —- | — = | ee ee 
BES aes ares: dc/ak, 473 lesa Ser? 


is the usual density of single-excitation states at the Fermi level. The 
integrals are over the Fermi surface. This result was first pointed out by 
Landau (1956), and its validity has been verified by the more fundamental 
calculations mentioned in § 2.6. Although in our nomenclature there are 
no interaction corrections to the specific heat, it should be noted that 
(4.1) is not equal to the specific heat in the Sommerfeld model because the 
band structure ¢«(k) contains exchange and correlation terms which affect 
the value of (4.2). What we mean is that the result (4.1) is given correctly 
by applying independent-electron formulae to the quasi-particle band 
structure e(k): we treat this band structure as the fundamental thing and 
the Sommerfeld model as a (convenient!) fiction of the mathematician’s 
imagination, so that we are only incidentally interested in the deviations 
from the Sommerfeld values. 


4.2. Paramagnetic Spin Susceptibility 

By this we mean the susceptibility due to turning some electron spins 
over, independent of spin-orbit or orbital contributions. Because we are 
altering the number of electrons with given spin, the total exchange energy 
is altered and this affects the spin susceptibility, as has been pointed 
out for a long time (Sampson and Seitz 1940, Pines 1955). Here we analyse 
the effect systematically in terms of our picture of excitations at the 
Fermi level. We consider the multiple excitation of fig. 12 with n electron- 
like excitations with up spin just about the Fermi level, and the same 
number of holes with down spin} just below e,. 

The electrons (holes) take up all states within an energy Se of the Fermi 
level, where 

n 

dA (ex) Q 
The one-particle contribution H, of (1.15) to the total energy is 


By =nb€ = 2n?|/N (e,)Q. 
Re een te 


+ We define a down-spin. hole-like excitation as one in which a down-spin 
electron has been removed, i.e. the system is left with a net up spin. 


be= 
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The interaction energy EH, we write as 
H,=3> >dn(k; * , k; @ )+3> dak, a kp ) 
+2 dalk,t , k, +) 


=-—vi,, 
where using (1.14) v can be cast in the form 
ier 02 e391 
MN (ep) N Ts (27r)6 
dS ds’ ‘ 
al jefok dejan AK TK T+ AKT ed DI. (43) 


By is negative, because A( + + ), consisting mostly of exchange energy, is 
negative and numerically larger than A(t, |). Thus the total energy, 
including magnetic, is 
ae an" 1—v) 
MN (ep) Q 
where pz is the Bohr magneton. By minimizing H with respect to n at 
constant H, we calculate the susceptibility y per unit volume as 


Xa Mmin/QH =WNV (ey) upr(1 = py. A ° : . (4.4) 


— 2npp 


Fig. 12 


Multiple excitation produced by magnetic field. 


The correction factor (1—v)~! is greater than unity because a negative H, 
makes it easier to turn spins over. Using the Bohm and Pines interaction 
for A, the factor is typically 1-2 to 1-5 (Pines 1955), but a more realistic 
interaction term would probably give values more like 1-1 to 1-2 as discussed 
in §5. 

Incidentally it should be noted that we give each excitation a magnetic 
moment of pz, not (e*/e)uy. The quasi-particle excitations were obtained 
in §1.2 by adding (subtracting) one electron to the system : it is irrelevant 
that only a fraction e*/e of this gets localized if one forms a wave 
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packet, the rest being uniformly distributed over the system, because 
the susceptibility is a bulk property and both the localizable and 
non-localizable parts contribute to the total magnetic moment. 


4.3. Anomalous Skin Effect 
The surface impedance of a metal in the extreme anomalous limit is 
governed by the limiting value as qo of the conductivity o(q) in a 
spatially varying electric field (Pippard 1960, eqns. (39) and (42)) 
6 (2) = 6 og EXP(— 192). 
More precisely we consider wavelengths 1/q small compared with the mean 
free path J, but large compared with the interatomic spacing so that the 
electric field does not see just the bare electron but the whole quasi-particle 
including all the correlation and exchange hole. 


Fig. 13 


> 


Excitation of Fermi surface in anomalous skin effect. 


The displacement of the Fermi surface at any point on the surface is 
(Pippard 1960, eqn. (15)) 
_ (1/h)re*(k)é, 


A= aa MS eink gil 5" 


For large qi (4.5) is large only over an angular range of order 6, 


around the positions P, P’, where v,(k)=0, Physically, only those electrons 
moving nearly parallel to the surface contribute much to the current. 
and this defines an effective zone PP’ (fig. 13) round the Fermi surface. 
We proceed to calculate the constants « and f in (1.20). Consider the 
total energy (1.15), and let us for a moment idealize (4.5) slightly as (fig. 13) 


A(k)=Ap; |v, /v| < 4, 
=0, — |v,/e]> 4. 
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The one-particle energy is approximately 


Ey (1/48) | 2Le(K)max — €p ]Aop,20 dk, 


= (1/429) [ (cfd, )y(Aetlat, 
=O(gl) 


The interaction energy is 


Bx (1/1679) | | i(k, k’)(Agp,20 dke,)(Ao’p,/20 dk’) 
= O(gl). 


Thus in the limit q/- oo, H, is negligible compared with H,, and it can easily 
be verified that this applies also with the correct displacement (3.5). 
We could now continue to evaluate (1.20). However, since the interaction 
energy H, is negligible, the correct answer must be given by the usual 
independent-particle analysis provided we give the quasi-particles the 
correct charge e*(k). The argument of Pippard can easily be modified 
in this way, and we obtain the surface impedance Z,, for diffuse surface 
reflection (Pippard 1960, eqns. (16) and (42)) 


B= Y=™ (1 +V/3i), 


S 


a (4.6) 
3— — | [e*(k)]*[p,| dx, 
=| [e(k)Peyldhy 
where p, is the radius of curvature at the effective point P in the section 
k,=const (fig. 13). Alternatively we can derive the same result using 
(1.20). 


4.4. Hlectrical Resistivity 


We have already derived in principle a formula for the d.c. resistivity 
R in §1.3.4 and it only remains to express the energy factor « and the 
current factor B in (1.20) as integrals over the Fermi surface. As in § 4.3, 
the acceleration of a quasi-particle is proportional to e,*, and the multiple 
excitation corresponds to a sideways shift of (e,*/e)A at each point of the 
Fermi surface in the direction of the electric field where A is now a constant. 
The total energy and current are then expressed as integrals over the Fermi 
surface. We shall not write down the explicit results, since the methods 
have already been illustrated in §§4.2 and 4.3 and the assumption of a 
constant relaxation time is so crude as to rob the results of any valid 
application. In §5 we shall discuss taking an arbitrary collision integral 
into account. 

However, it is interesting to analyse the case of free electrons further. 
By ‘ free electrons’ we mean electrons with coulomb interaction but without 
an external periodic potential. In this case the total force on the electron 


P.M.S. B 
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gas in a field & is Ne&, the acceleration of the centre of mass Veé/Nm, 
and the rate of increase of current 


2 _Neé x, 
accel 


dt ~ Nm 


At equilibrium this has to equal 


a a) 
dt collisions T 


and we obtain for the resistivity R 


peo 


This formula is therefore valid in the presence of coulomb interactions 
and does not depend on the kind of kinetic arguments found in most books. 
We shall now verify that it follows from (1.20) in our formalism. The 
electric field exerts the same force e& on each electron so that the whole 
system is accelerated without altering the relative motion of the electrons. 
Now the total energy / can be written, in quantum as in classical mechanics, 
E=E (motion of centre of mass)+H (motion relative to centre of mass). 


We therefore have 
E=P?/2Nm+ Eg, (4.8) 


where H, is the ground-state energy and P the total momentum. Also 


(§ 3.2) 
J=(e(m) Py” a A 
and from (1.16), (1.17), (1.20) we recover (4.7). 

Since (4.7) is also the result for the Sommerfeld model, does this mean 
the resistivity involves no interaction corrections? No it does not. By 
‘corrections’ in the context of this paper we mean corrections to the 
results obtained by applying the independent electron calculations to the 
band structure ¢(k) of the quasi-particles. This is the band structure which 
for instance is relevant to the electronic specific heat as shown in §4.1_ 
and includes exchange and correlation energy e,.(k); in our case 


e(k) = A2k2/2m4+ «(kj > 5k ae 440) 


Now in calculating the energy of the current carrying state, it is essential 
to include the interaction energy H, of (1.15) because it exactly cancels off 
the contribution of ¢,,(k) to the one particle energy H,, as can easily be 
shown in a general way. We obtain 


= (h?/2m)[ Sk? — Yk?) 
which is equivalent to (4.8). Similarly in calculating the current, we have 


v(k) A hk/m and e*(k) #e because of ¢,.(k), but taking the product e*(k)v(k) 
restores the simple result (4.9) as discussed in § 3.2. 
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With more complicated band structures, similar cancellations will occur 
so that it might well be sensible to analyse data in terms of the independent 
electron Eaotiel using a bogus band structure «’(k) which has the correct 


Fermi surface but has d¢’/0k adjusted in a way that corresponds to dropping 
e. from (4.10). 


4.5. Cyclotron Resonance 


In a metal, cyclotron resonance is usually observed with the magnetic 
field H parallel to the surface and the skin depth very small compared 
with the orbit radius (Pippard 1960). The quasi-particles are accelerated 
each time they come within the skin depth, and the resonant frequency 
is determined by the time taken to complete one revolution. The accelera- 
tion in k-space is 

e*(k) 


k = ae a 41 
k 5c v(k),AH ( ) 


Our reasons for including e* in this are given in §3.3. The motion of the 
representative point k is perpendicular to v(k) and to H, i.e. along an 
energy contour in a plane perpendicular to H. The time taken to traverse 
an element dk is hcdk/e*v, H where v, is the component of v perpendicular 
to H, and the periodic time of the orbit 


he [ dk 
Hero. 
The cyclotron frequency is therefore 


27H = 
= ai am CS 
ae a | bias k)v, (Kk | sal 


where the integral is taken round a section of the Fermi surface perpen- 
dicular to H. There are no interaction energy corrections to this, because 
the acceleration in the field is determined by the current j =e*v and there 
are no interaction corrections to the current as shown in §1.3.5. As it 
stands expressed in terms of the band structure, (4.12) involves e* rather 
than e, but for free electrons (with coulomb interactions) it reduces to the 
Sommerfeld value w,.=eH/mc on account of (3.25) as we might have 
deduced also by general arguments similar to those in § 4.4. 


4.6. De Haas—van Alphen Effect 


In a magnetic field, various properties including the diamagnetic 
susceptibility oscillate with constant period in 1/H. In the usual inde- 
pendent-electron model, the quantized electron orbits are given by 
(Pippard 1960, eqns. (11) and (25)) 


A =(n+ >)2reH |h, i parte S's we 4713) 
and the period of the oscillations by 
27 
= sO whey Cie a lie ace” Terk ee 


Here . is the area of the orbit in k-space, ,) the extreme cross section of 


P,M.S. I 
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the Fermi surface perpendicular to H, n and integer, and ¢ a phase factor. 
Unfortunately there appears to be no plausible way of modifying this 
analysis to take account of an effective charge: even assuming e* is a 
constant, we cannot with any conviction simply replace e by e* in (4.13), 
(4.14). . 

Let us therefore consider a different approach. The first question is: 
is the de Haas—van Alphen effect essentially a property of the whole 
ground state, or can we hope to discuss it by quantizing the excitations at 
the Fermi level? Although at first sight it would appear to be a ground- 
state property, two reasons suggest it may be correct to analyse it in 
terms of our excitations. Firstly in the independent electron model all 
the electrons are quantized according to (4.13), but for electrons well 
inside the Fermi surface this only amounts to a unitary transformation 
which does not affect the whole wave function, and just the quantization 
at the Fermi level is significant. Secondly, if we think of the temperature 
variation of the de Haas—van Alphen effect or the specific heat rather 
than the de Haas—van Alphen effect itself, it seems clear that the period 
is governed by the availability of quantized states for the thermally 
excited quasi—particles at the Fermi surface. 

We must therefore quantize the orbits of the quasi-particles. We have 
already derived in (4.12) the cyclotron frequency w, of revolution. Now 
from the correspondence principle, we expect the semi-classical absorption 
of energy in cyclotron resonance described in §4.5 to be equivalent to 
the jumping of quasi-particles to higher levels: ie. we expect fw, to 
equal the difference in energy between successive quantized levels, as is 
indeed the case in the independent electron theory (Pippard 1960, p. 193). 
This immediately implies that there should be some e* correction to the 
de Haas—van Alphen period, for if we quantize the orbits according to 
(4.13), the energy difference between successive levels as measured by 
e(k) is not hw, (4.12). We therefore look for a generalization of (4.13) 
which does give hw, (4.12) for the difference between levels. If we assume 
e*(k) to be independent of energy «(k) at each point on the Fermi surface 
but varying over the Fermi surface, we easily find 


l ( kadk 
which gives for the period 
27 [1 (kadk] 2 
A(1/H) = = sf | Ot a ara 
— : ext 


However, in general e* does vary with e, and there does not seem to be a 
general formula for the quantized orbits which when differentiated with 
respect to e gives (4.12); for such a formula would have to contain e* to 
get e* into (4.12) and thus would at the same time give terms in de*/de. 
We can only state that the period should involve some e*/e correction. 
Clearly the chain of reasoning is somewhat tenuous, and we must report 
that Luttinger (1960) profoundly disagrees with our result. His results 
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on the oscillations in the specific heat and the work of Idrimura and 
Tanaka (private communication) on the de Haas-—van Alphen effect, 
indicate that the period is correctly given by (4.14) without any e¢* correc- 
tion. It is not easy to see how this squares with the correspondence 
principle, as already discussed. 


4.7. Ultrasonic Attenuation 


The ultrasonic attenuation shows oscillations in a magnetic field when 
successive wave fronts of the ultrasonics measure the size of the quasi- 
particle orbits like a pair of calipers (Pippard 1960). In the case of constant 
e*, the sections of the Fermi surface have the same shape as the orbits 
in real space, but are larger by a factor e*H/h and rotated through 7/2, 
so that the ultrasonic oscillations enable the shape of the Fermi surface 
to be studied, but a factor e*/e is involved in the absolute size. If e* 
varies over the Fermi surface, the relation between the orbits in k-space 
and real space breaks down and no general relation seems possible. 

The ‘fork’, ‘tilt’ or “surf riding” effect (Reneker 1959, Spector 1960) 
on the other hand measures the Fermi velocity directly 

Us= v(Q)cos b 
were v, is the ultrasonic velocity, ¢ the critical angle between its direction 


of propagation and H, and v(Q) is the Fermi velocity at the point Q on 
the Fermi surface where the velocity is parallel to H. 


§ 5. Summary, DiscussIon AND COMPARISON WITH EXPERIMENT 


In §§1 and 2 we have described the low-energy excited states of an 
electron gas in terms of electron-like and hole-like excitations, where 
these quasi-particles have a well-defined band structure ¢(k) near a sharp 
Fermi surface given by «(k)=e«,. For a free electron gas this model seems 
to be onasound basis. In the presence of a periodic potential the existence 
of such excitations has not been proved yet, but it is difficult to believe 
that the same general picture does not apply because of the experimental 
evidence (§ 1.1) as well as physical reasonableness. An important property 
of the quasi-particles is that they carry a charge e*(k) which is not equal to 
the bare electronic charge e (§3). They also have a long lifetime which 
tends to infinity at the Fermi surface. Thus in discussing transport 
phenomena, we assume that a given transport process produces a multiple 
excitation of the system which can be represented by a certain distribution 
of quasi-particles, as shown for instance in fig. 1 for a uniform electric 
field. It has not been proved that the electric field acting on the whole 
system really produces a state like this. However, such states form semi- 
stationary states of the system with a long lifetime, so that if something 
different were produced by the electric field it would rapidly decay into 
the semi-stationary states. In any case, a multiple excitation of the form 
of fig. 1 would appear to be demanded by the way the quasi-particles are 
accelerated in the field, though it is a little obscure how the electric field 
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creates the elementary excitations in the first place. In fact the validity 
of our whole model needs verifying by more fundamental and rigorous 
calculations of particular properties. 

We have therefore a description of the electron system in terms of a gas 
of quasi-particles. Under almost all conditions the density of quasi- 
particles is low, of order &7'/e ~ 10~? times the electron density. For this 
reason they can be treated to a first approximation as independent. 
However, there is a screened coulomb repulsion energy between them, and 
the consequent interaction energy makes a significant contribution to any 
process involving energy balance or transfer, e.g. the rate of Joule heating 
in d.c. conduction (§ 1.3.4). The interaction also leads to a certain amount 
of scattering which can contribute to the electrical and thermal resis- 
tivities (§ 1.3.2). 

In order to focus attention on the role of the energy, electrical resistances 
are calculated from the rate of dissipation of energy as Joule heat. ‘The 
method of § 1.3.4 is limited by the assumption of a constant relaxation 
time 7. It is also physically oversimplified in that the primary mechanism 
destroying the current is usually elastic collisions which themselves do 
not dissipate the energy. The correct thing is to calculate the rate of 
entropy production by the primary collisions: in the steady state this 
entropy is passed around among the electrons before it finally appears in 
the lattice as heat. With this modification, our method becomes identical 
with the method of general transport theory based on a variational principle 
(Ziman 1960, p. 283). We see now that in order to take a general collision 
process into account, we only need a formula for the entropy of a given 
multiple-excitation including interaction energy to replace eqns. (7.8.4) 
to (7.8.6) of Ziman (1960). 

Something must be said now about the size of effects described. Un- 
fortunately experimental data are meagre and detailed calculations have 
only been made for the Bohm and Pines interaction (2.35). Consider 
first e*/e. The values calculated for a free electron gas with the Bohm and 
Pines interaction have already been given in table 2, and we note for example 
for r,=3 (corresponding approximately to copper) that e*/e=0-88. 
However, we expect that the Bohm and Pines interaction overestimates 
the effect considerably, and that the smoother (and more realistic) inter- 
actions shown in fig. 9 will give values of e*/e much nearer unity. The 
difference between hk,/m and v, comes from the rate of change of the 
exchange energy in ¢«(k) with k: an interaction which is a constant in 
k-space would give e*/e=1. To obtain a rough estimate of the effect, 
let us write for the mean value and the fluctuation of the interaction 


Vip =0-6 + 0-4, \ 
ewer estimate (k) =0:3 “3 O-1 


where the units and V(best estimate) are those of fig. 9. Since e*/e 
depends on the fluctuation of V(k), we expect V,)(k) to overestimate 
(e*/e)—1 by a factor of say 4, and a more realistic estimate is e*/e x 0-97 


(5.1) 
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As regards the sign of (e*/e)—1, ordinarily we expect the exchange 
energy to decrease as one goes to higher energy away from the Fermi 
surface, which gives e*/e<1. However, in real metals with complicated 
band structure effects, the possibility of e*/e>1 may not be completely 
ruled out. Experimentally the only measurements at present which 
allow us to determine e* appear to be those on the anomalous skin effect 
in copper and silver. The measurements (Pippard 1957, Morton 1960) 
were analysed using (4.6) with e* put equal to e. The apparent Fermi 
surface obtained in this way should therefore have a volume containing 
(e*/e)* instead of one electron per atom, if as a first approximation we 
assume e* is constant over the Fermi surface. Experimentally the volumes 
of the apparent Fermi surfaces were 1-00 + 12% (copper) and 0:95+9% 
(silver) electrons per atom. The results are therefore : 

For copper e*/e=1-00+ 0-04, 

For silver  e*/e=0-98 + 0-03. 
Thus the data are consistent with e*/ex 0-97 but do not give any proof 
that e*/e41. 

A comparison of (4.1) and (4.4) shows that a value for the correction 
factor (l1—v)-' may be obtained directly from measurements of the 
electronic specific heat and spin susceptibility : 


y $b? 
We use the data on sodium from Schumacher and Slichter (1956) for y 
and two recent measurements of y by Gaumer and Heer (1969) and Lien 
and Phillips (1960), giving 
(l1—v)+=1-20 (using Gaumer and Heer), 
=1:10 (using Lien and Phillips). 


The difference between these may be taken as an indication of the probable 
error, largely due perhaps to the partial phase transformation at low 
temperatures. The value given by the Bohm and Pines theory is 1-48: 
because of the double integration over the Fermi surface in (4.3), the 
value is more sensitive to the mean strength of the interaction than to its 
variation with k and we conclude from (5.1) that V,,, overestimates the 
effect by a factor of say 2, giving a more realistic estimate of (1—v)"' 1-24 
which is more nearly in agreement with experiment. 

Finally we turn to what is probably the greatest gap in the theory at 
present. Although the acceleration of the quasi-particles in an electric 
field appears to be correctly given—the theory is certainly self-consistent— 
uncertainties remain about the motion in a magnetic field, particularly as 
regards the de Haas—van Alphen effect. We have not been able to resolve 
theoretically whether or not the de Haas—van Alphen period should 
contain a correction factor involving e*/e. The experimental data are 
analysed using (4.14). Assuming (4.16) to be correct, this would give an 
apparent Fermi surface containing (e* je)? electrons per atom for a 
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monovalent metal. Experimentally Shoenberg (1960) finds the periods of 
the ‘belly’ orbits in copper to be consistently less than calculated for 
free electrons, giving the total volume of the apparent Fermi surface to 
be about 5% +14% too high. Allowing for the volume of the * ‘necks 

would only increase thi further. Thusif we assume there is no fundamental 
error in the calibration of the measurements, we would conclude firstly 
that the de Haas—van Alphen period does indeed contain a correction 
factor of e*/e and secondly that e*/ex0-97 + 0-01 for copper. 
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AE PEN? DE 


From (3.23) the number of electrons travelling in one quasi-particle is 
6N=Q lim > «*(k,)o(k;+q)B(k; k;t+q) . . . (Al) 


q>0) ki 
where the symbols are defined in §3.1. In this expression a(k;+q) can be 
chosen to be continuous and we can therefore put 
lim o(k;+ q) =a(k;), 
q>0 
and B(k;, k;+q) is a continuous function of k; as can easily be seen from 
its definition. Due to the small region in which a(k;,) is defined, we can take 
the value B(k;, k;+q)=B(ky, ky +q) and assume it is constant over the 
small volume in k-space where a(k;) 40. Using (3.7) for the normalization 
of the wave packet, (A 1) is transformed into 
1 ine . > A Re 
SN O lim (ky, ke aie. cee ee) 
q>U 
which is the limit we now evaluate. It is not incidentally equal to B(ko, ky) 
which is just the mean density N(+1)/Q. There is a genuine discontinuity 
in this function at q=0. To calculate the limit we study the commutator 


AP = p? H — Hp”, ‘ ; é F : ; (A 3) 
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where p°? is given by (3.5) and W by (2.1). Using the usual anticommuta- 


tion rele tions. of the operators a,*, a,, it follows after a long but straight- 
forward calculation 


= x {*.q+397}4,*a,, exp(iq.r). eo ed) 


We now use the definition (A3) and assume that k, is sufficiently near the 
Fermi surface that we can treat |k;) as an eigenvector of #: 


(ko, tA. IKo +4, t) =[e(k)+q) — e(Ky ) Pro, kO+q 
is exp (iq.r 24 k 
= oars) q. %.G+ 37} (Ko, tla,*a,, , glk +4, ¢). 
(45) 


Taking q#0 and using the Taylor expansion for «(ky +q) as well as the 
definition (3.11), (3.12) of B(ky, ky +q), (A5) reduces to 


Geo a feere 
(3,24) 0 mton 


h? } 
-a{= 3 > #4 kola,*4,,. glko +4) + 34B(Ky, kota} . (A6) 
and from this, using (3.2) and (A 2) 


v(k,)dNV 
= lim (h/m) > x kpla,*a,., 4|ko +4). ee ae AT) 


q>0 


It is lengthy but not difficult to prove that the function on the right-hand 
side is continuous for q=0 so that the limit is simply given by putting q=0. 
One way of proving the continuity is by expanding |k,) in the perturbation 
series 

[ky =A’|ky) + > A”(k,, k,)|k;, ko —k,;+k,; k,) +... . (A8) 


and similarly for |ky+q). Evaluating the right-hand side of (A7) we see 
that crossed terms (due to matrix elements between unperturbed states 
with different number of quasi-particles) are all at least proportional to q, 
and terms coming from matrix elements between states of equal number of 
quasi-particles can be expressed as 


(Lion) <ko|P lk; 3 Og) see te 8 (V9) 
where P is the momentum operator 
P=) xa,*a,.. ied Noh a LO) 
. 


Thus in the limit as q+0, we obtain from (A 7) 
V (le, )Ouv == (ale) (Kol PIKES), itt ae chars ee Rb) 


which proves that v is Sone to ras and by comparison with (3.18), 
(3.19) that 6N =e*/e. 
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Similarly to prove that the background gives no contribution to the 
current, we evaluate the expectation value of the current density operated 
for the wave packet (3.6) : 


I(r, t) = Ck, o, tJ (r)| ky; of) 2 ee ey 
jerry DY ua,*d,.q EXP (14-4). Le LS) 
md” a 
We obtain 
J(r, th= > ak (kja(ky) A 5 ej ot, Sa eee ae 
ki, kj 


he ? 
dict hea oe aa exp iq(r—v(k;)t | 
x > «¢ k;|a,,*a,., g|ki +4). ay St Shere Cad) 


We have proved before that the summation on the right-hand side of 
(A 15) is a continuous function of q. Consequently the coefficients of the 
Fourier series (A 14), (A 15) are continuous for the value q=0 and the 
current density packet represented by (A 14) is completely localized and 
similar to fig. 11(a). Therefore no contribution to the total current 
comes from the background and the spatial transference of charge is only 
made by the localized quasi-particle. 
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